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Array  Processing  for  Interference  Suppression  in  GPS 

Receivers 


Executive  Summary 

This  report  presents  research  work  performed  under  the  contract  no.  N000 14-04- 1-0607  with  the 
Office  of  Naval  Research  over  the  period  of  May  2004  to  September  2005.  It  includes  contributions  from 
Prof.  Moeness  Amin  (PI),  Dr.  Wei  Sun  (Postdoctoral  Fellow),  and  Mr.  Liyu  Liu  (graduate  Student). 
The  report  consists  of  two  chapters.  Each  chapter  has  its  own  abstract  and  introduction,  and  has  also  its 
own  equation  numbers,  figure  numbers,  and  references. 

In  Chapter  1,  we  consider  interference  suppression  and  multipath  mitigation  in  Global 
Navigation  Satellite  Systems  (GNSSs).  In  particular,  a  self-coherence  anti-jamming  scheme  is 
introduced  which  relies  on  the  unique  structure  of  the  coarse/acquisition  (C/A)  code  of  the 
satellite  signals.  Because  of  the  repetition  of  the  C/A-code  within  each  navigation  symbol,  the 
satellite  signals  exhibit  strong  self-coherence  between  chip-rate  samples  separated  by  integer 
multiples  of  the  spreading  gain.  The  proposed  scheme  utilizes  this  inherent  self-coherence 
property  to  excise  interferes  that  have  different  temporal  structures  from  that  of  the  satellite 
signals.  Using  a  multiantenna  navigation  receiver,  the  proposed  approach  obtains  the  optimal  set 
of  beamforming  coefficients  by  maximizing  the  cross  correlation  between  the  output  signal  and  a 
reference  signal,  which  is  generated  from  the  received  data.  It  is  demonstrated  that  the  proposed 
scheme  can  provide  high  gains  toward  all  satellites  in  the  field  of  view,  while  suppressing  strong 
interferes.  By  imposing  constraints  on  the  beamformer,  the  proposed  method  is  also  capable  of 
mitigating  multipath  that  enters  the  receiver  from  or  near  the  horizon.  No  knowledge  of  either  the 
transmitted  navigation  symbols  or  the  satellite  positions  is  required. 
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In  Chapter  2,  we  study  the  performance  of  the  delay  lock  loops  (DLL)  in  Global  Positioning 
Systems  (GPS)  receivers  in  the  presence  of  non-Gaussian  noise.  GPS  receivers  find  growing 
applications  in  indoor  and  outdoor  communication  environments  including  urban  and  rural  areas. 
Interference  and  noise  sources  for  GPS  receivers  may  assume  Gaussian  or  non-Gaussian 
distributions,  specifically  when  operating  inside  buildings.  The  GPS  receiver  performance  under 
Gaussian  additive  noise  has  been  extensively  studied.  Non-Gaussian  noise  may  equally 
contaminate  the  GPS  satellite  signals  and  compromise  the  receiver  delay  lock  loops  (DLL), 
producing  significant  tracking  errors.  These  sources  include  impulsive  noise,  ultra-wideband 
(UWB)  signals,  and  impulse  and  noise  radar  signals  for  target  tracking  and  through-the-wall 
imaging  applications.  This  chapter  considers  non-Gaussian  noise  and  examines  its  effect  on  GPS 
correlation  and  discriminator  outputs  for  the  commercial  GPS  using  Coarse  Acquisition  (C/A) 
code.  The  correlator  noise  output  components  are  produced  from  the  correlation  between  the 
noise  sequence  and  the  early,  late,  and  punctual  reference  C/A  code.  Due  to  the  long  time¬ 
averaging  performed  in  the  GPS  correlation  loops,  these  components  assume  Gaussian 
distributions.  The  discriminator  tracking  error  variance  is  derived,  incorporating  the  effect  of 
noise,  the  front-end  precorrelation  filter,  and  the  sampling  rate. 
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Chapter  1 


A  Novel  Interference  Suppression  Scheme  for  Global 
Navigation  Satellite  Systems  Using  Antenna  Array 

1.  Introduction 

Satellite  navigation  is  a  tool  to  determine  position,  velocity,  and  precise  time  world  wide.  A  navigation 
receiving  device  determines  its  three  dimensional  position  plus  time  by  measuring  the  signal  traveling 
time  from  the  satellite  to  the  receiver  (the  so  called  pseudorange  due  to  the  clock  offset  at  the  receiver) 
[1],  The  generic  name  of  the  satellite  navigation  systems  is  Global  Navigation  Satellite  System  (GNSS). 
Currently,  there  are  two  operative  navigation  systems,  one  is  Global  Positioning  System  (GPS)  of  the 
United  States  and  the  other  one  is  Global  Navigation  Satellite  System  (GLONASS)  of  Russia  [2], 

Both  GPS  and  GLONASS  employ  direct  sequence  spread  spectrum  (DS  SS)  signaling.  The  GPS 
satellites  transmit  signals  at  two  L  band  frequencies:  LI  =  1:57542GHz  and  L2  =  1:2276GHz.  Each 
satellite  broadcasts  two  different  pseudorandom  (PRN)  codes,  a  coarse/acquisition  (C/A)  code  and  a 
precision  (P)  code,  using  code  division  multiple  access  (CDMA)  technique.  TheLlcarrier  transmits  both 
the  C/A  code  and  the  P  code,  whereas  the  L2  carrier  only  transmits  the  P  code  [1],  The  GPS  C/A  code  is  a 
Gold  code  with  a  chip  rate  of  1.023  Mchips/sec  (or  code  period  1023)  and  repeats  every  millisecond,  and 
the  P  code,  usually  encrypted  for  military  use,  has  a  chip  rate  10.23  Mchips/sec  and  repeats  about  every 
week.  Similar  to  GPS,  GLONASS  also  has  two  DS  SS  components.  However,  frequency  division 
multiple  access  (FDMA)  technique  is  used  in  GLONASS,  where  each  satellite  transmits  on  a  different 
center  frequency.  The  C/A  code  of  GLONASS  has  a  length  of  51 1  chips  at  a  chip  rate  of  51 1  kHz,  and  it 
is  the  same  for  every  GLONASS  satellite.  The  C/A  code  repeats  10  times  within  each  navigation  symbol 
which  has  a  rate  of  100  bps.  Another  component  of  GLONASS  has  10  times  the  chip  rate  (5.11  MHz)  of 
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the  C/A  code  and  uses  a  longer  PRN  code.  In  this  chapter,  we  consider  only  the  signals  induced  by  the 
C/A  code. 

Despite  the  ever  increasing  civilian  applications,  the  main  drawback  of  the  satellite  navigation 
systems  remains  to  be  its  high  sensibility  to  interference  and  multipath  [3],  [4],  [5],  [6],  [7],  [8],  which  are 
the  two  main  sources  of  errors  in  range  and  position  estimations.  The  effect  of  interference  on  the  GNSS 
receiver  is  to  reduce  the  signal  to  noise  ratio  (SNR)  of  the  navigation  signal  such  that  the  receiver  is 
unable  to  obtain  measurements  from  the  satellite,  thus  losing  its  ability  to  navigate  [1],  Jammers  reported 
to  impact  the  GPS  receivers  are  wideband  noise,  CW,  pulsed  noise,  pulsed  CW,  swept  tone  (chirped), 
frequency  hopping,  and  spoofers.  Each  type  of  jammers  has  advantages  and  drawbacks  in  terms  of 
complexity  and  effectiveness  [9],  The  spread  spectrum  (SS)  scheme,  which  underlines  the  GNSS  signal 
structure,  provides  a  certain  degree  of  protection  against  interference.  However,  when  the  interferer  power 
becomes  much  stronger  than  the  signal  power,  the  spreading  gain  alone  is  insufficient  to  yield  any 
meaningful  information.  For  example,  for  the  GPS  C/A  signal,  the  receiver  is  vulnerable  to  strong 
interferers  whose  power  exceeds  the  approximately  30  dB  gain  (101og101023~30dB)  offered  via  the 

spreading/despreading  process.  It  is  thus  desirable  that  the  GNSS  receivers  operate  efficiently  in  the 
presence  of  strong  interference,  whether  it  is  intentional  or  unintentional. 

Interference  suppression  in  SS  communication  systems  has  been  an  active  research  topic  for  many 
years  and  a  number  of  techniques  have  been  developed  (see,  e.g.,  [4],  [10],  [11],  [12],  and  references 
therein).  In  satellite  navigation,  interference  can  be  combated  in  the  time,  space,  or  frequency  domain,  or 
in  a  domain  of  joint  variables,  e.g.,  time  frequency  [13],  [14]  or  space  time  [15],  [16],  Multiple  antenna 
receivers  allow  the  implementations  of  spatial  nulling  and  beamsteering  based  on  adaptive  beamforming 
and  high  resolution  direction  finding  methods.  These  methods  are  considered  to  be  effective  tools  for  anti 
jam  GPS  [9]. 

Conventional  antenna  arrays,  which  are  only  based  on  spatial  processing,  are  among  the  simplest,  and 
yet  effective,  techniques  for  narrowband  interference  suppression.  Such  techniques,  however,  is 
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inadequate  for  broadband  jammers  (such  as  spoofer)  cancellation  or  in  the  presence  of  multipath.  In  these 
cases,  the  temporal  degree  of  freedom  is  required.  Space  time  processing  provides  the  receiver  with 
spatial  and  temporal  selectivity.  The  spatial  selectivity  allows  the  discrimination  between  the  navigation 
and  interference  signals  based  on  their  respective  direction  of  arrivals  (DOA’s)  [17],  [18],  [19],  [20],  The 
temporal  selectivity  is  used  for  broadband  interference  and  multipath  cancellation.  Generally,  the  criteria 
for  determining  the  optimal  array  weights  include  maximum  signal  to  interference  plus  noise  ratio 
(MSINR),  minimum  mean  square  error  (MMSE),  and  minimum  output  power  (MOP)  [16].  The  MSINR 
approach  seeks  the  array  weight  vector  by  maximizing  the  receiver  output  SINR.  The  MMSE  method 
chooses  the  weight  vector  such  that  the  mean  square  difference  between  the  array  output  and  the  desired 
temporal  signal  is  minimized.  Since  the  navigation  signal  power  is  well  below  the  noise  floor  at  the 
receiver,  minimizing  the  output  power  while  attempting  to  preserve  the  navigation  signal  is  the  goal  of  the 
MOP  based  scheme.  While  these  methods  are  widely  used  in  interference  suppression  in  satellite 
navigation  systems,  one  obvious  drawback  is  that  they  all  require  some  kind  of  a  priori  knowledge  of  the 
problem  parameter  values.  For  example,  satellite  locations  are  needed  in  order  to  calculate  the  signal 
power  for  the  MSINR  and  MMSE  methods.  In  [21],  spatial  and  temporal  processing  techniques  are 
applied  to  remove  GPS  like  broadband  jammers  and  recover  the  navigation  information.  The  assumptions 
made  in  [21]  are  that  the  chip  and  bit  synchronizations  are  achieved,  implying  that  pseudorange 
measurements  are  obtained.  However,  the  assumption  of  satellite  positions  or  acquisitions  is  difficult  to 
enforce  under  persistent  jamming,  or  during  the  initial  satellite  searching  stage  when  any  synchronization 
is  yet  to  be  established. 

In  addition  to  interference,  GPS  pseudorange  and  carrier  phase  measurements  also  suffer  from  a 
variety  of  systematic  biases,  including  satellite  orbit  prediction  error  and  clock  drift,  ionospheric  and 
tropospheric  delay,  GPS  receiver  clock  offset,  and  signal  multipath  [22].  The  satellite  orbit  and  timing, 
ionospheric,  and  troposhperic  errors  can  be  removed  by  differencing  techniques  or  significantly  reduced 
by  modeling  [1].  The  receiver  clock  offset  can  also  be  corrected  by  differencing  but  is  often  solved  for  as 
an  unknown  in  the  position  solution.  Multipath,  on  the  other  hand,  is  normally  uncorrelated  between 
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antenna  locations.  As  a  result,  differencing  will  not  cancel  the  errors  caused  by  multipath.  Also,  modeling 
multipath  for  each  antenna  location  is  difficult  and  impractical  [23].  To  combat  signal  multipath,  many 
techniques  have  been  proposed.  Among  them,  narrow  correlator  is  one  of  the  most  widely  used 
approaches  that  improves  the  C/A  code  tracking  performance  by  reducing  the  space  between  the  early  and 
late  correlator  [24].  Other  multipath  mitigation  techniques  include  multipath  elimination  delay  lock  loop 
(MEDLL)  [25]  and  multipath  estimation  technology  (MET)  [22],  etc. 

This  chapter  proposes  a  new  interference  suppression  technique  for  GNSS  using  spatial  processing, 
but  incorporating  the  known  temporal  structure  of  the  C/A  signal.  A  careful  examination  of  the  existing 
interference  cancellation  techniques  reveals  that,  though  efficient  in  most  situations,  they  do  not  fully  take 
advantage  of  the  unique  C/A  signal  structure,  namely  the  replication  of  the  C/A  code.  Due  to  the 
repetition  of  the  spreading  code,  the  GNSS  C/A  signal  exhibits  strong  self  coherence  between  chip 
samples  that  are  separated  by  integer  multiples  of  the  spreading  gain.  Utilizing  this  feature,  an  anti 
jamming  technique  is  developed  to  suppress  a  large  class  of  narrowband  and  broadband  interferers.  It  also 
has  the  capability  of  mitigating  multipath,  resulting  in  improved  accuracy  in  pseudorange  measurements. 
The  proposed  technique  allows  the  civilian  C/A  code  tracking  and  acquisition  operations  in  the  presence 
of  strong  interference,  specically  at  "cold  start”,  where  there  is  no  prior  information  on  satellite  angular 
positions  or  ranges  [26].  In  military  applications,  the  encrypted  P  code  is  used  instead  of  the  C/A  code. 
However,  due  to  the  short  duration  of  the  P  code  chip  (10  times  shorter  than  the  C/A  code  chip),  the 
synchronization  in  P  code  is  usually  difficult  to  achieve  using  the  early  late  correlator,  and  assistance 
from  the  C/A  code  is  needed  [24].  With  the  receiver  introduced  in  this  chapter,  initial  synchronization  in 
the  P  code  can  be  established  by  first  processing  the  interference  suppressed  C/A  code.  In  essence,  the 
self  coherence  based  anti  jamming  approach  is  a  blind  technique,  which  does  not  require  the  knowledge 
of  the  navigation  data  or  satellite  locations  to  perform  interference  suppression.  This  makes  it  most 
applicable  in  the  initial  satellite  searching  phase  when  such  information  is  unavailable,  or  in  a  prolonged 
jamming  environment  where  the  formerly  obtained  satellite  positions  are  no  longer  reliable. 
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The  rest  of  the  chapter  is  organized  as  follows.  In  Section  2,  we  briefly  review  the  concept  of  spectral 
self  coherence  restoral.  An  anti  jamming  GNSS  scheme  is  presented  in  Section  3  and  two  receivers  based 
on  such  scheme  are  developed.  In  Section  4,  we  discuss  various  important  issues  related  to  the 
performance  of  the  proposed  GPS  receiver.  Numerical  results  are  presented  in  Section  5  to  demonstrate 
the  performance  of  the  proposed  receiver.  Finally,  Section  6  concludes  this  section. 


2.  Overview  of  SCORE  Algorithms 


The  proposed  anti  jamming  GNSS  technique  builds  on  the  basic  concept  of  the  self  coherence  restoral 
algorithm  proposed  in  [27].  A  signal  s(f)  is  referred  to  as  spectrally  self  coherent  at  frequency  separation  /? 
if  the  correlation  between  the  signal  and  its  frequency  shifted  version  is  nonzero  for  some  lag  x,  i.e.,  if 
[27],  [28] 


(0.1) 


where  (*)*  denotes  the  complex  conjugate  and  represents  the  infinite  time  averaging 

operation.  p{sf >  (t)  is  the  self-coherence  function  and  Rss  (0)  and  (t)  represent  the  average  power 
and  cyclic  autocorrelation  function  of  5(f),  respectively.  For  an  M-e lement  vector  waveform  x(r),  the 
cyclic  autocorrelation  matrix  R^(t)  is  defined  as 


Rif  (T)  □  (x(t)xH  (1 t  -  T)e-j2Kp' )_  (0.2) 

where  (~)H  stands  for  the  complex  conjugate  transpose.  Complex  wide  sense  cyclostationary  waveforms 
exhibit  spectral  self  coherence  at  discrete  multiples  of  the  time  periodicities  of  the  waveform  statistics 
[27].  The  signal  waveforms  that  possess  the  self  coherence  feature  include  most  communication  signals, 
such  as  PCM  signals  and  BPSK  signals  [28]. 
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The  spectral  self  coherence  restoral  (SCORE)  beamforming  techniques  have  been  shown  to  blindly 
extract  the  desired  signal  in  the  presence  of  unknown  noise  and  interference  [27].  The  SCORE  algorithms 
seek  the  beamformer  weight  vector  that  maximizes  a  measure  of  the  cyclic  feature  of  the  beamformer 
output.  For  example,  in  the  presence  of  interference,  the  received  signal  is  given  by  x  (/)  =  as(t )  +  v  (*), 
where  a  is  the  signal  amplitude  and  v  (t)  is  the  interference,  which  is  assumed  to  be  independent  of  s(f).  If 
£(/)  is  spectrally  self  coherent  at  a  frequency  shift  /?,  then  the  cyclic  autocorrelation  of  x(t)  can  be 
expressed  as  [27] 

(r)  =  H2  Rf  (r)  +  R(nfi)  (t)  =  \af  Rjf>  (t)  (0.3) 

Equation  (1.3)  shows  that  the  shift  in  frequency  completely  decorrelates  the  interference  component  in 
x(0,  given  that  v(f)  is  not  spectrally  self-coherent  at  the  frequency  separation  /?. 

There  are  several  different  versions  of  the  SCORE  algorithm,  of  which  the  least-squares  (LS)  SCORE 
is  the  simplest.  The  LS-SCORE  algorithm  determines  the  array  weight  vector  by  minimizing  the 
difference  between  the  array  output  and  a  reference  signal,  which  is  obtained  by  processing  the  delayed 
and  frequency-shifted  version  of  the  received  signal.  Other  SCORE  algorithms  include  the  cross-SCORE 
algorithm,  which  determines  the  beamformer  by  strengthening  the  cross-correlation  between  the  output  of 
the  array  and  a  reference  signal,  and  the  auto-SCORE  algorithm,  which  maximizes  the  spectral  self¬ 
coherence  strength  at  the  output  of  a  linear  combiner  [27].  The  self-coherence  anti-jamming  scheme 
proposed  in  this  section  is  based  on  the  cross-SCORE  algorithm. 

3.  Proposed  Anti-Jamming  GNSS  Scheme 

Before  introducing  the  proposed  anti-jamming  receiver,  we  first  examine  the  temporal  structure  of  the 
navigation  signal,  as  the  receiver  is  developed  by  exploiting  the  repetitive  feature  of  the  C/A  signal. 
Figure  1  depicts  the  structure  of  the  received  noise-free  navigation  signal,  where  the  BPSK  modulated 
navigation  symbols  (simply  referred  to  as  “symbol”  thereafter)  are  spread  by  a  PRN  code  with  spreading 
gain  of  P  (P  =  1023  for  GPS  and  P  =  511  for  GLONASS)  and  chip-rate  sampled.  The  code  sequence 
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(denoted  as  “spreading  block”  in  Figure  1)  is  repeated  L  times  ( L  =  20  for  GPS  and  L  =  10  for 
GLONASS)  within  each  symbol.  Two  blocks  of  data  are  formed  at  the  receiver:  a  data  block,  which 
spans  N  consecutive  samples,  and  a  reference  block  with  the  same  number  of  samples  as  the  date  block. 

The  distance  between  the  respective  samples  in  the  data  and  reference  blocks  is  set  equal  to  jP  chips, 

/ 

where  1  <  j  <  L  .  Obviously,  due  to  the  repetition  of  the  spreading  code,  the  nth  sample  in  the  data  block 
will  have  the  same  value  as  the  corresponding  nth  sample  in  the  reference  block,  providing  that  the  two 
samples  belong  to  the  same  symbol. 

From  the  temporal  structure  of  the  C/A  signal,  we  observe  the  inherent  self-coherence  between 
samples  in  the  data  block  and  the  reference  block,  due  to  the  repetition  of  the  spreading  code.  Based  on 
this  observation,  a  novel  anti-jamming  technique  is  developed  in  this  section,  and  Figure  2  shows  a  block 
diagram  of  the  proposed  receiver  with  this  technique.  In  Figure  2,  an  M-element  array  is  deployed.  There 
are  two  beamformers  in  the  receiver:  a  main  beamformer  w,  processing  samples  in  the  data  block,  and  an 
auxiliary  beamformer  f,  handling  data  from  the  reference  block.  An  error  signal  e(t)  is  formed  as  the 
difference  between  the  beamformer  output  z{t)  and  a  reference  signal  d(t ),  which  is  the  output  of  f  .  For 
the  proposed  scheme,  the  weight  vectors  w  and  f  are  updated  according  to  the  cross-SCORE  algorithm. 
The  signal  reaching  the  GNSS  receiver  is  the  aggregate  of  the  satellite  navigation  signals,  their  respective 
multipaths,  additive  white  Gaussian  noise  (AWGN),  and  broadband/narrowband  interferers.  Thus,  after 
carrier  synchronization,  the  signal  received  at  the  receiver  can  be  expressed  as 

x(«)  =  £  sq(nTs  -  Tq )aqeJ%  +£  Bkuk  (n)dk  +  v(n) 

<7=0  *=  1 

K  (0.4) 

=  s0(nTs  - TQ)a0eM‘  +YJsq(nTs  ~ Tq > V*  +  Z +  v(«) 

g= 1  *=1 

where  Ts  is  the  Nyquist  sampling  interval,  Q  is  the  number  of  multipath  components,  with  subscript  0 
designated  to  the  direct-path  signal.  Due  to  the  weak  cross-correlation  of  the  C/A-codes,  only  one  satellite 
is  considered  in  Equation  (1.4).  In  the  above  equation,  sq(ri),  T  ,  and  <j)q  are  the  signal  sample,  time- 
delay,  and  phase-shift  of  the  qth  multipath  component,  respectively,  K  is  the  number  of  interferers,  uk{n) 
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is  the  waveform  of  the  ki h  interferer  with  amplitude  Bk.  The  vectors  and  d*  are,  respectively,  M*  1 
spatial  signatures  of  the  qth  satellite  multipath  and  the  Mi  interferer,  and  v(/i)  consists  of  noise  samples. 

Let  s(n )  □  s0(nTs  —T0)a0ej^°  denote  the  data  vector  across  the  array  due  to  the  direct-path  signal.  Then, 
Equation  (1.4)  can  be  rewritten  as 

x(n)  =  s(n) + s(n )  4-  u(rc)  +  \{n)  (0.5) 

Q  ..  K 

where  s(/z)  □  ^sq(nTs  —  Tq)aqe}^  and  u(n)  0  Bkuk (n)dk  .  Assuming  a  direct-path  to  the  satellite  at 

9=1  *=i 

direction  /y  and  a  uniform  linear  array,  we  can  express  vector  ao  in  the  specific  format  of  a  steering  vector 
as 

a0=a(#)Cl[l,  ej2*fcT,  •••,  jr  (06) 

where  fc  is  the  carrier  frequency,  T  —  A/csin#  is  the  interelement  path  delay  of  the  source  in  the 
direction  of  0,  c  is  the  propagation  speed  of  the  waveform,  and  A  is  the  sensor  spacing.  According  to  the 
formulation  of  the  data  and  reference  blocks,  the  counterpart  of  x(n)  in  the  reference  block  within  the 
same  symbol  can  be  written  as 

X(n  -  jP)  =YjSg  (nTS  ~Tq-  jP>qe^q  +Y.BkUk(n~  jp)dk  +  V(«  -  jP) 

9=0  *=1 

=2>, (nTs  -  tq )a/^  +  £ Bkuk ( n -  jP)dk  +  v(n  -  jP)  (0.7) 

q= 0  k=\ 

=  s(n)+s(n)+u(n-  jP)  +  \(n  —  jP) 
where  we  have  assumed  that,  when  considered  within  the  same  symbol, 

sq{nTs-Tq)  =  sq{nTs-tq-  jP),  q  =  0,  •••,  Q  (0.8) 

Compared  to  the  general  case  of  self-coherence,  there  is  no  frequency  difference  between  the  signal 
samples  in  the  data  and  reference  blocks,  i.e.,  the  frequency  shift  /?  =  0. 
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From  Figure  2,  the  beamformer  output  and  the  reference  signal  are  given  by 

d(n)  □  iHx(n  -  jP) ,  respectively.  We  define  the  following  covariances: 

z{n)  □  w Hx(n)  and 

Rzd  0  E{z(n)dH(nj\  =  wHE{x(n)xH(n-jP)}( 

(0.9) 

R a  □  E{z(n)zH  (n)}  =  wH E{x(n)xH  (n)}w 

(0.10) 

Rdd  D  E{d(n)dH(n)}  =  fHE{x(n-jP)xH(n-jP)}f 

(0.11) 

Under  the  assumption  that  the  navigation  signal,  interference,  and  noise  are  independent,  then 

R  „UE  {x(n)xH  («)}  =  £{x(n  -  jP)xH  (n  -  jP )}  =  RS+RU  +  RV  (0.12) 


The  three  terms  on  the  right-hand  side  of  Equation  (0. 12)  denote,  respectively,  the  covariance  matrices 
of  the  C/A  signal,  including  both  the  direct-path  and  multipath  signals,  interference,  and  noise: 

Rs  0£{[s(n)  +  s(n)][s(n)  +  s(n)f},  R„  □  £{u(n)u" («)} ,  and  Rv  □  £{v(n)v"(n)}  .  If  the 

navigation  signal  is  the  only  data  component  which  correlates  with  its  delayed  version,  then  the  cross¬ 
correlation  matrix  between  the  corresponding  data  vectors  in  the  data  and  reference  blocks  simplifies  to 
[cf.  Equation  (0.3)] 

R^  0  E{x(n)xH(n-jP)}  =  R,  (0.13) 

3.1.  Cross-SCORE  Algorithm  Based  Receiver 

We  first  consider  the  receiver  design  by  directly  applying  the  cross-SCORE  algorithm.  For  the  proposed 
GNSS  anti-jamming  scheme,  there  are  two  beamformers  w  and  f  to  be  determined.  With  d(n)  serving  as 

the  reference  signal,  we  define  e(n)  □  z(n)—d(n)  as  the  difference  between  the  receiver’s  output  and 
the  reference  signal.  The  relationship  between  w  and  f  can  be  established  in  the  least-squares  (LS)  sense. 
For  a  fixed  beamformer  w,  the  LS  solution  of  f  is  given  by  fLS  =  R^r^ ,  where 

=  E{x(n-yP)z//(n)]  =  R(pHw.Similarly,  if/ is  fixed,  then  =R^R^)Wf . 
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According  to  the  cross-SCORE  algorithm,  the  beamformers  w  and  f  are  obtained  by  maximizing  the 
cross-correlation  between  z(n )  and  d(n): 


C(w,f)  □ 


K\ 

2 

w"RLP)f 

2 

dd 

>"R„w] 

(0.14) 


Substituting  f  and  w  in  the  above  equation  by  and  Wls,  respectively,  we  have 


C(w,fu) 


w^R^R-’R^w 


wwRv 


vv 


(0.15) 


fgR<?R;X?f 

fwR,  f 


(0.16) 


The  weight  vectors  w  and  f  that  maximize  C(w,  fLs)  and  C(Wls>  f),  respectively,  are  readily  shown  to  be 
the  eigenvectors  corresponding  to  the  largest  eigenvalues  of  the  generalized  eigenvalue  problems: 

R„»  =  2R«R;XP,w  (0.17) 

Rjf  =  (0.18) 


where  X  and  k  are  the  eigenvalues. 

It  is  observed  from  Equations  (0.17)  and  (1.18)  that,  for  the  proposed  receiver,  the  main  beamformer 
w,  which  generates  the  receiver  outputs,  is  equivalent  to  the  auxiliary  beamformer  f  that  provides  the 
reference  signal.  This  equivalence,  however,  is  not  surprising  because  of  the  unique  structure  of  the  C/A 
signal.  From  Equations  (0.4)  and  (1.7),  we  note  that  the  self-coherence  of  the  navigation  signal  is  due  to 
the  time  lag  between  the  two  samples  which  do  not  encounter  any  frequency  shift  after  the  frequency 
demodulation.  Therefore,  x(n)  and  x(n-jP)  have  the  same  correlation  function,  given  by  Equation  (1.11). 
For  the  general  case  of  self-coherence,  on  the  other  hand,  the  signal  auto-correlation  function  does  not 
necessarily  equal  to  the  auto-correlation  function  of  the  frequency-shifted,  time-lagged  version  of  the 
original  signal.  As  a  result,  the  cross-SCORE  algorithm  will  not  produce  two  identical  beamformers  [27], 
as  it  does  for  the  proposed  receiver. 
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Simulations  presented  in  Section  V  show  that  the  cross-SCORE  based  receiver  performs  fairly  well  in 
a  jamming  environment.  It  is  capable  of  producing  high  gains  for  satellites  currently  in  the  field  of  view, 
while  suppressing  strong  jammers.  However,  such  a  receiver  provides  no  measures  against  multipath, 
which  is  one  of  the  dominant  error  sources  in  navigation.  To  address  the  multipath  issue,  we  modify  the 
receiver  design,  as  discussed  in  the  next  subsection. 

3.2.  Modified  Cross-SCORE  Algorithm  Based  Receiver 

Since  the  GNSS  signal  multipath  shares  the  same  structure  as  the  direct-path  signal,  it  is  expected  that 
both  the  C/A  signal  and  the  undesired  multipath  components  will  appear  at  the  receiver’s  output 
undistorted. 

It  is  known  that  the  satellites  typically  lie  above  the  horizon,  whereas  the  multipath  is  often  generated 
from  local  scatters  near  the  horizon  [16].  To  equip  the  receiver  with  means  to  combat  multipath,  while 
maintaining  the  self-coherent  approach,  we  introduce  constraints  on  f  such  that  the  reference  signal  d(n ) 
does  not  contain  reflections  from  near  the  horizon.  To  do  so,  we  define  equally  spaced  directions 

d=l,  •*,£>,  covering  some  solid  angle  Q.  near  the  horizon.  Let  B  □  [b(ft)  •••  b (Yd)\  be  the  M  xD 


matrix  consisting  of  steering  vectors  defined  as  in  Equation  (6).  To  mitigate  multipath  in  the  range  O,  we 
require  =  0.  Then,  the  cost  function  in  Equation  (1.16)  is  rewritten  as 

f hR  f 

fop,  =  ar§ max  "  >  subject  to  Bwf  =  0  (0. 19) 

f  1  Kxr‘ 


where  .  The  solution  of  the  above  equation  is  obtained  as  follows.  Let  r  =  rank(B) 


<  min(M,D)  be  the  rank  of  the  B  matrix.  Performing  the  singular  value  decomposition  (SVD)  [29]  of  B 
yields 


U"B"V 


A  0"| 

Lo  dj 


(0.20) 


where  U  and  V  are  two  unitary  matrices  with  dimension  MxM and  D x  D,  respectively,  and 
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A  =  diag{<7,,  <j2,  •••  crr} 


(0.21) 


where  (J,  >  <72  >  •  •  •  >  <Tr  are  eigenvalues  of  B  arranged  in  a  decreasing  order.  Let  A  be  formed  from  the 
last  M-r  columns  of  U.  Thus,  A  spans  the  null  space  of  Bw,  i.e.,  BWA  =  0.  Let  a  be  a  (M-r)  x  1  vector  such 


f  =  Aa  (0.22) 

Using  vector  a,  the  constrained  maximization  problem  in  Equation  (1.19)  is  transformed  to  an 
unconstrained  one: 

awAwR  Aa 

a  =  arg  max  — - — 5 — 25 —  (0.23) 

&  «  a^R^Aa 

The  above  generalized  eigendecomposition  problem  can  be  solved  using  Cholesky  decomposition. 
Particularly,  since  R„  is  positive  definite,  A^R^A  is  also  positive  definite.  Then  the  Cholesky 
decomposition  of  AHR„A  is  A7/R„A  =  GGW,  where  G  is  a  (M-r)  x  (M-r)  matrix  with  full  rank  [29]  and, 

thus,  invertible.  Let  a  =  G_wp  .  Then 

«YR  Aa  PWG  ‘A^R  AG  Hp  p"G'‘A"R AG'Hp 

XX  _  1  XX  ■  ■  XX  *  /a  Ayi  \ 

a^A^R^Aa  ~  p^G  ^GG^G  ^p  pwp 

Accordingly,  the  maximization  problem  becomes 


a^R^Aa  p/,G,AwRxcAG/,p 

max— — - — — =  max- - rH2 - -  =  maxp  G  A  R  AG  P  (0.25) 

«  a*AHR  Aa  *  * 


xx 


under  the  standard  constraint  lipil  =  1,  where  INI  is  the  vector  2-norm  [29].  Hence,  P  is  given  by  the 
eigenvector  associated  with  the  maximum  eigenvalue  of  G_1  AwR;ccAG_H  .  And,  finally. 


fc  =AGWP 


The  beamformer  w  is  derived  as 
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according  to  the  LS  relation  between  w  and  f. 

4.  Covariance  Matrix  Estimations 

In  practice,  the  covariance  matrices  Rxf  and  R(J/  are  unknown  and  have  to  be  replaced  by  their  sample 
estimates.  Define  the  M  xN  data  and  reference  matrices  as  XN  □  [x(n),  •••,  x(n  —  (N -{))]  and 

XNte{  D[x(n  —  jP),  x(n  —  (N  —  1)  —  jP)] ,  where  N  is  the  block  length  and  N  <P.  The  sample 

covariance  matrices  are  then  given  by 

k,  <°-28> 

R«=Tx„Xl,  (0.29) 

And,  the  beamformers  w  and  f  are  calculated  correspondingly. 

It  is  noted  from  Equation  (1.16)  that  the  covariance  matrices  R^and  determine  the 

performance  of  the  proposed  receiver.  In  practical  implementations,  the  data  and  reference  blocks  and 
X/tof  are  used  to  estimate  R^  and  R^} ,  and  subsequently  provide  the  weight  vector  w,  which  is  then 
applied  to  process  signal  samples  in  the  data  block.  The  key  assumption  made  for  the  proposed  GPS 
receiver  in  Section  HI  is  that  both  the  data  and  reference  samples,  x(n)  and  x(n-  jP),  1  <j  <  20,  belong  to 
the  same  navigation  symbol.  However,  since  the  data  samples  used  for  covariance  matrix  estimations  are 
selected  randomly,  and  interference  suppression  is  performed  prior  to  any  symbol  synchronization 
process,  there  is  no  guarantee  that  the  data  and  reference  samples  belong  to  the  same  symbol.  Questions 
arise  as  how  will  the  receiver  perform  when  the  above  assumption  fails,  i.e.,  the  data  and  reference 
samples  lie  in  two  adjacent  symbols? 
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To  answer  the  above  question,  we  relax  the  condition  imposed  in  Section  3,  and  develop  the  general 
expression  of  the  covariance  matrices  between  the  data  and  reference  samples,  x(n)  and  x(n-yP).  Define 
the  following  events: 

A :  x(n)  &  x(n  -  jP)  are  within  the  same  symbol, 

A  2 :  x(ri)  &  x(n  -  jP)  are  in  two  adjacent  symbols, 

(0.30) 

A  2| :  x(n)  &x(n-  jP)  are  in  two  symbols  with  the  same  sign, 

A  22 :  x(n)  &x(n-  jP)  are  in  two  symbols  with  different  signs 

With  random  selection  of  time  n,  and  using  the  repetitive  property  of  the  C/A-code,  it  is 
straightforward  to  show  that  the  corresponding  probabilities  of  the  above  events  are 

Pr{Ai}  =  -J-  = 1  ~|r .  Pr(A2 }  =  ^T>  Pr{A2i}  =  ^>and  Pr{A22}  =  ^,  respectively,  where  T 

=  20 P  is  the  total  number  of  samples  in  one  symbol.  The  exact  expression  of  the  cross-correlation 
function  R^  can  be  written  in  terms  of  the  above  probabilities  and  conditional  expectations  as 

R^=£{x(«)xlf(»-yP)|A1}ft:{AI}  +  £{x(»)xir(n-7P)|A2}Pr{A2}' 

=  Ejxfnjx^^-y'PjlAjPrjA,}  (0.31) 

+E  {x(/z)xw  ( n  -  ;'P)|  A  21 }  Pr  {A  21 }  +  E  {x(n)xH  (n  -  y'P)|  A  22}  Pr  {A  22} 

Since  £{x(n)x"(n- jP)|a,}  =  P{x(n)xH(n- yP)|A21}  =  R,  and 

E{x(n)xH  ( n  -  ;P)|A22}  =  -R  s ,  then 

RL"  =R,(Pr{A,}+Pr{A2l})-R,Pr{A22}=fl-^lR,  (0.32) 

Equation  (1.32)  shows  that  the  covariance  matrix  R^  depends  on  the  distance  between  the  data  and 

reference  samples  jP.  The  maximum  value  of  is  achieved  when  j  =  1,  representing  the  closest 
possible  data  and  reference  blocks. 
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In  practice,  however,  sample  estimates,  obtained  from  Equations  (1.27)  and  (1.28)  using  the  data  and 
reference  blocks  XN  and  XWref,  replace  the  exact  values  in  Equation  (1.32).  It  can  be  readily  shown  that  if 
X;V  and  X,Vl,,f  are  jP  samples  apart,  1  <j  <  20,  the  probability  of  the  two  blocks  belonging  to  the  same 

jP+N 

symbol  or,  equivalently,  in  two  adjacent  symbols  with  the  same  sign,  is  1 - .  On  the  other  hand, 


jP-N 

the  probability  that  XN  and  XWref  are  in  two  adjacent  symbols  with  opposite  signs  is  — — — .  Using  the 


above  probabilities,  the  expected  values  of  Rxt  and  are  derived  in  Appendix  A  as: 


R  +R  +R 


(0.33) 


which  show  the  same  dependency  on  jP  as  in  Equation  (1.32)  and  that  R^  and  R^?  are  unbiased 
estimates  of  R^  and  R^J ,  respectively. 

The  above  covariance  matrix  estimations  use  only  one  data  block  and  its  replicated  reference  block. 
To  fully  take  advantage  of  the  repetitive  feature  of  the  C/A-code,  multiple  data/reference  blocks  can  be 
used  in  the  time-averaging.  Particularly,  using  G  data  and  reference  blocks.  Equations  (1.27)  and  (1.28), 
respectively,  become 

*«o=4lX(!r>x!!<*>/w  <°-35> 

G  g=l 

RS=^Ex»(s)JOs>/JV  <»■»> 

G  g=1 

In  the  case  when  one  of  the  data  blocks  (and,  respectively,  a  reference  block)  is  split  between  two 
adjacent  symbols  with  opposite  signs,  a  maximum  of  only  two  of  the  G  terms  in  the  above  equation  may 
suffer  from  symbol  transition,  whereas  the  rest  of  the  terms  will  be  coherently  combined.  Appendix  A  and 
B  derive  the  mean  and  variance  of  the  above  estimations,  showing  the  value  of  using  a  higher  value  of  G. 
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5.  Numerical  Results 


In  this  section,  we  evaluate  the  performance  of  the  proposed  self-coherence  anti-jamming  receiver  using 
the  GPS  C/A  signals. 

A  uniform  linear  array  (ULA)  consisting  of  M  =  7  sensors  with  half-wavelength  spacing  is  used  in 
simulations  with  one  satellite  and  no  multipath.  We  set  M  =  11  for  simulations  with  multiple  satellites  or 
multipath.  The  GPS  navigation  symbols  are  in  the  BPSK  format  and  spread  by  C/A-codes  (Gold  codes) 
with  processing  gain  of  P  =  1023.  We  select  the  first  satellite  C/A-code  for  concept  demonstration.  At  the 
receiver,  chip-rate  sampling  is  performed  and  N  =  800  samples  are  collected  in  both  the  data  and 
reference  blocks  for  covariance  matrix  estimations.  The  signal-to-noise  ratio  (SNR)  and  signal-to- 
interference-plus-noise  ratio  (SINR)  are  defined,  respectively,  as  SNR  =  101og10l/Pvand 

SINR  =  101og10  1  /(/)  +  PV) ,  all  in  dB,  where  unit  signal  power  is  assumed  for  simplicity,  Pv  is  the  noise 

power,  and  P,  is  the  interference  power.  The  jammer-to-signal  ratio  (JSR)  is  defined  as  JSR  =  101og10  Pt 

dB.  Interferers  used  in  the  simulations  are  generated  as  broadband  binary  signals  having  the  same  rates  as 
the  C/A-codes,  but  with  a  different  structure  than  that  of  the  C/A  signals. 

5.1.  Antenna  Beam  Pattern  without  Interference 

We  first  consider  the  scenarios  in  which  no  interference  presents  at  the  receiver.  SNR  is  -30  dB. 
Covariance  matrices  are  estimated  using  one  data  block  and  one  reference  block  taken  within  the  same 
symbol.  The  performance  of  the  cross-SCORE  based  receiver  is  shown  in  Figure  3,  where  the  antenna 

pattern  is  formed  towards  the  satellite  located  at  6  =  30° . 

We  recall  the  discussions  in  Section  4  which  suggest  that  better  performance  can  be  expected  when 
multiple  data  and  reference  blocks  are  used  to  estimate  the  covariance  matrices.  In  Appendix  B,  the 

variances  of  the  sample  estimates  of  are  calculated  and  it  shows  that  using  multiple  data/reference 
blocks  can  indeed  reduce  the  estimation  variance.  We  now  demonstrate  experimentally  the  effect  of 
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multi-block  estimation  on  the  receiver  performance.  Particularly,  a  very  special  situation  is  created  where 
one  of  the  data  blocks  is  evenly  split  between  two  symbols  having  opposite  signs.  SNR  is  set  at  -40  dB.  If 
G  =  2  and  the  split  data  block  happens  to  be  the  second  one,  the  receiver  fails  to  provide  any  substantial 

gain  for  the  satellite  located  at  30°,  as  shown  in  Figure  4(a).  This  is  because  that  elements  in  the  time¬ 
averaging  of  R(i'’).  given  by  Equation  (0.35)  cancel  each  other,  resulting  in  significantly  weak  cross¬ 
correlation  between  z{n)  and  d(ri)  [cf.  Equation  (0.14)].  If,  on  the  other  hand,  G  >  2  data  and  reference 
blocks  are  involved  in  the  estimation,  the  split  of  one  block  will  not  have  such  a  dramatic  impact  on  the 
receiver  performance  as  in  the  G  =  2  case,  as  only  two  among  G  blocks  are  affected  due  to  the  split.  It  is 
clear  from  Figure  4(b)  that  a  beam  is  generated  towards  the  satellite  with  G  =  7  despite  the  split. 
Generally,  to  avoid  the  performance  degradation,  it  is  recommended  that  odd  number  of  samples  (N) 
should  be  chosen  for  each  data/reference  block  and  2  <G<  D. 

It  is  known  that  in  satellite  navigation,  at  least  four  satellites  are  needed  simultaneously  in  the  field  of 
view  in  order  to  calculate  the  receiver’s  three-dimensional  position  and  time.  Since  the  proposed  receiver 
relies  on  the  special  structure  of  the  C/A  signals  to  suppress  interference  and  all  satellite  emitted  C/A 
signals  share  the  same  repetitive  feature,  it  is  expected  that  the  receiver  will  pass  the  signals  from  all 

satellites  with  high  gains.  In  the  simulation,  the  satellites  are  located  at  dx  =10°,  0,  =30°,  03  =50°, 

and  04  =  70° ,  with  SNR  =  -30  dB.  As  shown  in  Figure  5,  four  clear  beams  are  generated  towards  the  four 
satellites. 

A  point  worth  mentioning  is  that  the  receiver  presented  in  this  section  is  able  to  suppress  interference 
for  all  satellites  at  once.  In  GNSS,  since  the  satellite  spreading  codes  are  known  at  the  receiver,  it  may  be 
intuitive  to  consider  using  the  spreading  code  as  the  reference  instead  of  generating  one  from  the  received 
data.  Even  though  the  locally  generated  spreading  code  is  noise  and  interference  free,  it  cannot  serve  as 
the  reference  signal  in  the  proposed  receiver  because  1)  the  alignment  of  the  incoming  signal  and  the  local 
reference  code  may  not  be  established  during  the  interference  removal  stage,  hence  there  is  no  guarantee 
that  the  data  block  and  reference  blocks  are  separated  by  integer  multiples  of  P  chips;  2)  it  is  not  possible 
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to  use  one  specific  satellite’s  spreading  code  as  the  reference  signal  to  remove  interference  for  all 
satellites.  Therefore,  interference  suppression  must  occur  in  a  serial  manner.  The  disadvantage  is  obvious 
as  compared  to  the  simultaneous  interference  removal  the  proposed  receiver  offers. 

5.2.  Interference  Suppression 

We  next  investigate  the  receiver’s  interference  suppression  capability  by  comparing  it  with  the  MMSE 
receiver  of  [16].  The  MMSE  receiver  determines  the  weight  vector  by  minimizing  the  mean  square 
difference  between  the  array  output  and  the  desired  signal.  The  latter  approach,  however,  requires  the 
knowledge  of  the  satellite  direction.  This  condition  is  eliminated  in  the  proposed  scheme.  If  the  jammers 
have  explicit  bearings,  we  can  generate  the  received  signals  according  to  Equation  (1.4),  but  replacing  the 
spatial  signature  d*  by  the  respective  steering  vector  defined  in  Equation  (1.6).  The  direction  of  the 

satellite  is  20° ,  while  two  jammers  are  located  at  40°  and  60° .  The  weight  vector  in  the  MMSE  method 
is  obtained  by  using  the  exact  transmitted  navigation  signal.  Figure  6  clearly  shows  that  deep  nulls  are 
placed  at  the  jammer  locations,  whereas  high  gains  are  generated  towards  the  direction  of  the  satellite  in 
both  schemes.  The  advantage  of  the  proposed  receiver  is  that  neither  prior  synchronization  nor  known 
satellite  location  is  required. 

5.3.  Multipath  Effects 

The  purpose  of  the  simulations  performed  in  this  subsection  is  to  demonstrate  the  difference  between  the 
receiver  that  solely  relies  on  the  cross-SCORE  algorithm  and  the  receiver  with  additional  constraint  in  the 
presence  of  multipath.  As  discussed  in  Section  3.1,  the  cross-SCORE  based  receiver  is  unable  to  mitigate 
the  signal  multipath,  though  very  efficient  in  suppressing  interference,  as  shown  by  simulations  presented 
so  far.  This  drawback  has  motivated  a  constrained  receiver  design  and  resulted  in  the  modified  cross- 
SCORE  algorithm  based  receiver  in  Section  3.2. 
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We  assume  that  multipath  reaches  the  receiver  from  the  15-degree  range  (Q  as  defined  in  Section  3.2) 
above  the  horizon.  We  divide  Q  into  seven  2-degree  spaced  angles  and  form  the  corresponding  matrix  B. 
The  power  of  the  multipath  component  is  one  fifth  of  the  direct-path  signal  power. 

We  first  consider  the  case  when  there  is  no  interference.  The  direct-path  signal  is  incident  on  the 

array  with  50°  angle.  One  multipath  component  (half-chip  relative  delay  and  half  the  direct-path  signal 

power)  arrives  from  the  9°  direction.  Using  the  cross-SCORE  based  receiver,  both  the  direct-path  signal 
and  the  multipath  component  receive  high  gain  at  the  receiver  output  [Figure  7(a)].  If,  instead,  the 
modified  cross-SCORE  based  receiver  is  employed,  the  multipath  contribution  is  significantly  reduced 
from  the  output  of  the  receiver,  which  is  evident  from  Figure  7(b). 

In  the  next  case,  a  jammer  enters  into  the  system  from  30° .  Figure  8  shows  how  the  two  different 
receivers  respond  in  this  environment.  We  note  from  Figure  8(a)  that  both  receivers  can  successfully 
place  deep  null  at  the  jammer  location.  However,  the  two  receivers’  responses  to  multipath  are  just 
opposite.  While  the  cross-SCORE  based  receiver  generates  a  beam  towards  the  multipath  component,  the 
modified  cross-SCORE  algorithm  based  receiver  creates  a  null  at  the  same  direction. 

The  multipath  mitigation  performance  of  the  proposed  receiver  is  also  evaluated  by  feeding  the  output 
of  the  receiver  to  a  conventional  early-late  delay  lock  loop  (DLL)  [24],  We  consider  the  discriminator 
functions  of  the  receiver  outputs  without  multipath  mitigation  and  the  outputs  with  the  modified  cross- 
SCORE  algorithm.  We  compare  the  results  with  the  case  where  there  is  no  multipath.  The  early-late 
spacing  is  set  to  be  half  of  the  C/A  chip  interval.  The  simulation  results  are  depicted  in  Figure  8(b),  which 
clearly  shows  that,  without  any  multipath  mitigation  process,  the  zero-crossing  point  of  the  discriminator 
function  drifts  away  from  the  origin,  indicating  the  pseudorange  measurement  error  [30].  If,  on  the  other 
hand,  the  modified  cross-SCORE  receiver  is  used  first  to  mitigate  multipath  contributions,  the  zero¬ 
crossing  point  of  the  corresponding  discriminator  function  almost  overlaps  with  the  zero-crossing  point 
obtained  using  the  direct-path  only  signal,  suggesting  that  the  proposed  technique  can  significantly  reduce 
the  multipath  effect  on  pseudorange  measurement. 
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These  simulations  prove  that  both  receivers  have  the  capability  of  canceling  strong  jammers. 
However,  for  multipath  mitigation,  only  the  modified  cross-SCORE  algorithm  based  receiver  can  reject 
multipath  coming  from  near  the  horizon. 

5.4.  Synchronization  Process 

In  satellite  navigation,  the  receiver  is  ultimately  evaluated  by  its  ability  to  provide  accurate  pseudorange 
measurements.  This  is  achieved  by  establishing  synchronization  between  the  receiver  and  the  satellite, 
which  is  decided  based  on  the  cross-correlation  between  the  beamformer  outputs  and  a  locally  generated 
spreading  sequence  [31].  When  the  phase  of  the  receiver  replica  code  matches  that  of  the  code  sequence 
emitted  from  the  satellite,  there  is  a  maximum  correlation.  The  high-gain  beams  towards  the  satellites 
provided  in  the  previous  examples  should  be  examined  in  the  context  of  their  effects  on  the  post¬ 
processing  pseudorange  calculations. 

In  the  simulation,  the  satellite  is  located  at  20°  and  the  two  jammers  are  at  40°  and  60° .  The  figure 
of  merit  is  the  cross-correlation  between  the  receiver  output  and  the  Gold  code  sequence: 

37) 

where  c  denotes  the  P  x  1  receiver  Gold  code,  z  is  a  P  x  1  vector  with  elements  given  by  z(n )  =  wHx 

and  w  is  the  beamformer  coefficient  vector  discussed  in  Section  3.1.  The  normalized  cross-correlation 
with  the  respective  antenna  beam  pattern  for  SNR  =  -25  dB  and  JSR  =  30  dB  and  50  dB  are  shown  in 
Figures  9  and  10,  respectively.  Also  shown  in  these  figures  are  the  normalized  cross-correlations  obtained 
before  the  jammers  are  removed.  It  is  observed  from  Figure  9(b)  that  synchronization  can  be  achieved  in 
the  presence  of  interference  when  JSR  =  30  dB.  Figure  10  shows  that  the  proposed  receiver  can 
effectively  cancel  directional  jammers  and  achieve  synchronization  even  when  the  JSR  is  as  high  as  50 
dB  [Figure  10(c)].  Without  interference  suppression,  however,  synchronization  fails  as  shown  in  Figure 
10(b). 
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5.5.  Circular  Array 


In  addition  to  the  uniform  linear  array,  we  also  implemented  the  proposed  receiver  with  a  uniform  circular 
array  (UCA),  whose  configuration  is  shown  in  Figure  1 1(a).  Let  ( 6 \  0)  denote  the  elevation  angle  and  the 
azimuth  angle  of  the  satellite.  Then,  the  steering  vector  of  the  satellite  for  the  M-element  UCA  is  given  by 


a(0,0)  = 


.llir  •  a  a 
j~ — sintfcos^ 
y  X 


.Inr  .  _  ( ,  _  M-iY] 


U 


(0.38) 


where  r  is  the  radius  of  the  circular  array  and  X  is  the  wavelength.  The  steering  vector  of  the  jammer  has 
the  same  form  as  a(#  0)  given  above.  In  the  simulation,  the  satellite  signal  reaches  the  array  from 

( 10° ,  20° ),  whereas  a  jammer  is  located  at  (  60° ,  40° ).  The  beam  pattern  is  shown  in  Figure  1 1(b)  for  r  = 
X,  and  SNR  =  -30  dB  and  JSR  =  30  dB.  It  is  observed  from  Figure  1 1(b)  that  the  receiver  has  the  ability  to 
reject  jammers  from  arbitrary  directions. 


6.  Conclusions 

In  this  chapter  ,  we  addressed  the  issue  of  interference  suppression  in  global  satellite  navigation  system. 
Specifically,  the  unique  structure  of  the  GNSS  C/A  signal  is  exploited.  Due  to  the  repetition  of  the  C/A- 
code  within  each  navigation  symbol,  strong  self-coherence  is  observed  between  chip-rate  sampled  signals. 
It  is  shown  that  the  use  of  this  self-coherence  feature  allows  the  development  of  an  anti-jamming  GNSS 
receiving  scheme  which  is  built  on  the  cross-SCORE  algorithm.  The  proposed  scheme  incorporates 
multiple  data  and  reference  blocks,  separated  by  integer  multiples  of  the  spreading  code  length,  to 
generate  the  array  weight  vectors.  Its  performance  is  analyzed  in  view  of  its  dependency  on  the  number 
of  blocks  and  the  number  of  samples  in  each  block.  Two  receivers  are  constructed  based  on  the  proposed 
scheme.  One  directly  applies  the  cross-SCORE  algorithm  which  seeks  the  optimal  beamformers  by 
maximizing  the  cross-correlation  between  the  receiver  output  and  a  reference  signal,  derived  from  the 
receiver  signal.  The  other  one  applies  constraints  on  the  beamformer  such  that  it  can  also  reject  multipath 
arriving  from  and  near  the  horizon.  Numerical  results  have  shown  that  the  proposed  scheme  is  capable  of 
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suppressing  strong  wide  class  of  narrowband  and  broadband  interferes  while  preserving  signals  and  no  a 
priori  knowledge  of  either  the  transmitted  signals  or  the  satellite  locations  is  required. 


Appendix  A  Mean  Calculation 

To  simplify  the  derivation,  we  rewrite  the  received  signal  vector  as 

\{n)  =  s(n)a(0)  +  v(n)  (A.  l) 

where  we  consider  only  the  direct-path  signal  with  an  explicit  direction  0of  the  satellite.  The  vector  v(n) 
contains  samples  of  interference  and  noise,  with  zero  mean  and  variance  erf .  Both  components  of  v(n) 
are  assumed  to  be  independent  of  the  GPS  signal.  Accordingly, 
Rs  =  E{s(n)a(0)aH  (0)sH  (n)}  =  a(0)aH(0)  □  Rfl,  where  it  is  assumed 

E{s(n)sH  (n)j  =  E  {|.s(n)|2j  =  1 . 

From  Section  3,  the  estimates  of  the  covariance  matrices  R  h  and  Rff’’  are  obtained  using  the  data 
and  reference  blocks  XN  and  XNref  as  RXi  =  XWX"  /n  and  R^)=  XnX%k{/n  ,  respectively.  Taking 

A 

the  expected  value  R ^  yields 

R„  =  £{«„}  =  E{x(n  -  i)xH(n  -  <)}  =  R,  +  R. 

**  1=0 

The  expected  value  of  R^)  is  obtained  as  follows.  Define  the  following  events: 

Bx :  XN  &  XWref  are  within  the  same  symbol, 

B 2  :XN  or  XWref  is  split  between  two  adjacent  symbols 
&  Af,<  N  samples  are  in  the  first  symbol, 

B3 :  the  entire  XN  and  XNk{  are  in  two  adjacent  symbols. 


(A.2) 


(A.3) 
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The  corresponding  probabilities  of  the  above  events  are:  Pr{ff,}  =  l  —  ,  Pr {B^  =  -—-f  and 

jP-N 

Pr{#3}  =  - —  — ,  respectively.  In  addition,  we  define  the  following  events  regarding  the  two  adjacent 


symbols: 


CJ :  the  two  adjacent  symbols  have  the  same  sign, 
C2 :  the  two  adjacent  symbols  have  opposite  signs. 


(A.4) 


The  expected  value  of  is  calculated  as 


RL”  =  £{R‘r} = jrt.  «{*(»- Ox"  (n-i-jP)} = +  r«, + R 

**  1=1 


(P) 

xrlSi 


(A.5) 


where 


1  N 


□-E£{x(n-/)xH(n-/-yP)|fi1}Pr{B1}  =  (  1- 

i-\ 


jP  +  N 
T  J 


R 


(A.6) 


Because  the  GPS  symbols  are  equi-probable,  then  the  occurrence  of  Bi  implies 


p(P) 


□  E {x(“ “ *): x" (n-i-  /P)| B3 } Pr {S3}  =  0 


(A.  7) 


On  the  other  hand,  when  event  B0  occurs, 


i  N 

U—^E{x(n-i)xH  (n-i-  jP)\B2]?r{B2}=~Rs  (A.8) 


N . 


From  Equations  (A.5  -  A.8),  the  expected  value  of  R^}  is  given  by 


R<?  = 


1- 


jP+N 


N 


R,+— -R.  =  1-—  R 


jP 


(A.9) 


which  is  exactly  the  one  shown  in  Equation  (1.32). 

When  using  G  data  and  reference  blocks,  the  estimate  of  the  covariance  matrix  R^  is  given  by 
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(A.  10) 


Cr  g=1 


^  /  p\  ^  /  p\  ^ 

To  calculate  Ru(i ,  the  expected  value  of  Rxt(i ,  we  define  the  following  events: 

F, :  the  first  data  block  or  the  last  reference  block  is  split, 

F|:  one  of  the  other  G  - 1  data/reference  blocks  is  split 

&  G,<  G  data  blocks  (respectively,  G,  - 1  reference  blocks)  are  in  one  symbol,  (A.i  l) 
F3 :  the  data  blocks  and  reference  blocks  are  within  the  same  symbol, 

F4  :  no  split  block  and  the  data  and  reference  blocks  are  in  two  adjacent  symbols. 

The  corresponding  probabilities  are  Pr{F,}  =  -^-,  Pr{F2}  =  ~ ,  Pr{F3}  =  1  —  —  ,  and 

Pr  {F4}  =  .  We  maintain  that 


RSS1F,=^i£{x„(«)XL(*)k,}Pr{F,}/V=|-^R,+|-R 


RS,f,=7ri;F{x»(«)XL,(^)|F2}Pr{Fj/«  = 
U  «=1 


(G-l)iVG-2R  |  (G-l)iV R 
IT  G  s  2T  s 


Further, 


£{x„  (OX^^F,}  Pr{F,}/  N  =  )r 


£{X«  (#>x^r  («>|  F-4}  Pr{F- J  /  Af  = 

^  ^=1 


2  T  G 


Finally,  is  given  by 


i>(p)  _»(P)  ,pw  .£</>)  .p(?)  _  i  Jr  n 

1Vju:G  “  IKxxG\Fl  "r  1VaaGIF2  ^  rVxtGlF3  ~r  1VxcGIF4  “  1  ^ 

V  1 
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which  is  equivalent  to  the  expected  value  given  in  Equation  (1.32). 


Appendix  B  Variance  Calculation 


The  variance  of  R^  is  given  by  [32] 


var[R«]  =  £{R«Rr}-£J{Rlr} 


In  Appendix  A,  we  have  shown  that  the  covariance  estimates  are  unbiased.  In  what  follows,  we 


concentrate  on  evaluating  E  j  . 


Using  one  data  and  reference  block,  we  let  vPllfii  □  j  Pr  {/?,},  where  events  B, , 


le  [1,3],  are  defined  in  Equation  (A.3).  Then, 


0  £{R"RLP'"}  =  +V  + (A.18) 

When  the  data  and  reference  blocks  are  within  the  same  symbol  or,  equivalently,  in  two  adjacent  symbols 
with  the  same  sign  (i.e.,  event  CJ ),  we  have 


tf|Ci}=772-jE;^SSx(w-Ox"(n-/-j/>)x(w-/-j/>)xW(«-/)|cq |  (A.  19) 


Substituting  x(n)  from  Equation  (1.39)  and  after  some  straightforward  calculations,  we  have 


(„„ui  M( N+^)v  .  m+al) 


£{ft«Rr|c}= 


R..0R 


where  we  have  used  aw  (0)a(0)  =  M  and  E |vw  (n)v(n)j  =  M<rl .  When  event  Bl  occurs, 


Similarly, 
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(A.  22) 


and  in  case  of  event  ,  we  have 


Finally,  the  variance  of  R^}  is 


var[RL'>]  =  £{RL,>Rr}-£2{^’} 


.jicw+oJ)  ,£<1±0„  Jj.zTr, 


For  G  data  and  reference  blocks  we  can  similarly  define,  using  the  events  in  Equation  (A.l  1), 

*0  D  «{RS*S"} = + v  + 


where  1„:nE{^R™  F;jPr{F;i},  Ze  [1,4].  Following  the  same  procedure  we  adopted  in 


calculating  the  expected  value  of  R^} ,  it  can  be  readily  shown  that 


V  - 

*■  n\  v 


r  g2  T'j 

(G— l)iV(G-l)2  [  (G-pl- 


2T  'J 


GQP-V)  (G-l)2  |  GQP-W)1ft 


IT  'J 


from  which  we  obtain 


g{«Sig}  =  fl-2— -f2^-TM(W  +  ff-)R.+-^(1+<T-)R,|1 

1  J  v  T  GT  J|_  AT  AT  ‘J 


Finally,  the  variance  is 
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(A.  28) 


var[R«]  =  £{RrRr}-£J{6»} 

J \-2iL+2JL\m±tiK +m±^  RJ_f  ,_z 

V  T  GTJ  N  “  N  vll  V  T 


\2 


R 


The  above  equation  clearly  shows  that  the  larger  the  number  of  data  and  reference  blocks  used  in  the 
time-averaging,  the  smaller  the  variance. 
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Figure  1.  Noise-free  C/A  signal  structure  and  data  and  reference  blocks  formation. 
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Figure  2.  Structure  of  the  proposed  anti-jamming  scheme. 
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Figure  3.  Beam  pattern  generated  by  the  cross-SCORE  based  receiver  \ 
and  one  reference  block  taking  within  the  same  symbol. 


Figure  4.  Beam  pattern  generated  by  the  cross-SCORE  based  receiver  with  multiple  data  and  reference 
blocks  and  SNR  =  -40  dB.  (a)  G  =  2;  (b)  G  =  7. 
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Figure  5.  Beam  pattern  generated  by  the  cross-SCORE  based  receiver  with  four  satellites,  SNR  =-30  dB, 
and  G  =  7. 


Figure  6.  Antenna  gains  of  the  proposed  scheme  and  the  MMSE  scheme  with  SINR  =  -33  dB,  JSR  =  30 
dB,  and  G  =  3. 
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Figure  7.  In  the  presence  of  multipath  with  SNR  =  -30  dB  and  G  =  7.  (a)  Cross-SCORE  based  receiver; 


(b)  Modified  cross-SCORE  based  receiver. 


Figure  8.  (a)  Comparison  between  the  cross-SCORE  based  receiver  and  the  modified  cross-SCORE  based 
receiver  SINR  =  -  33  dB,  JSR  =  30  dB,  and  G  =  1\  (b)  Comparision  of  the  discriminator  functions  of  the 
early-late  delay  lock  loop. 
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Figure  9.  Synchronization  of  the  cross-SCORE  based  receiver  with  SNR  =  -25  dB,  JSR  =  30  dB,  and  G 
1.  (a)  Beam  pattern;  (b)  Normalized  cross-correlation  before  jammer  removal;  (c)  Normalized  cross¬ 
correlation  after  jammer  removal. 
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Figure  10.  Synchronization  of  the  cross-SCORE  based  receiver  with  SNR  =  -25  dB,  JSR  =  50  dB,  and  G 
=  7.  (a)  Beam  pattern;  (b)  Normalized  cross-correlation  before  jammer  removal;  (c)  Normalized  cross¬ 
correlation  after  jammer  removal. 


40 


Satellite 
(10,  20) 


Figure  11.  Performance  of  the  receiver  with  a  circular  array,  (a)  Circular  array  configuration;  (b)  Beam 
pattern. 
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Chapter  2 

Performance  Analysis  of  GPS  Receivers  in  Non- 
Gaussian  Noise  Incorporating  Precorrelation  Filtering 

and  Sampling  Rate 

1.  Introduction 

The  GPS  receiver  is  subject  to  noise  and  jamming,  in  addition  to  multipaths.  The  effect  of  multipath  on 
the  GPS  receiver  performance  has  been  extensively  analyzed  [1-5],  whereas  several  anti-jam  GPS 
techniques  have  been  proposed  for  suppressing  narrowband,  wideband,  noise-like  and  chirp-like  jammers 
[6-8],  This  chapter  focuses  on  effect  of  various  noise  sources  on  GPS  receiver  performance. 

With  the  growing  needs  for  indoor  GPS  and  demands  of  the  satellite-based  navigation  in  manned 
and  autonomous  ground,  aerial,  and  surface  vehicles,  the  GPS  receivers  are  required  to  operate  in  close 
proximity  to  various  noise  sources  of  Gaussian  and  non-Gaussian  characteristics.  These  sources 
compromise  the  receiver  performance  by  impeding  signal  acquisition  and  increasing  the  bit  error  rates 
[9].  Sources  of  non-Gaussian  noise  may  include  motor  ignition  noise,  which  is  generated  by  spark  plugs 
used  in  internal  combustion  engines  [10],  impulsive  and  noise  radars  for  imaging  [11],  and  ultra- 
wideband  (UWB)  emitters  [12,13].  The  latter  is  becoming  the  choice  of  signals  in  imaging  and  wideband 
communication  platforms  [14-18].  The  UWB  sources  and  noise  radar  have  bandwidths  that  can  extend 
over  HF,  VHF,  UHF,  and  beyond,  including  the  GPS  band.  The  transmission  of  repeated  short  pulses  in 
UWB  systems  contributes  to  their  possible  impulsiveness  nature.  Noise  radars  [11],  are  seeing  increased 
applications  in  indoor  and  outdoor  target  detection  and  tracking.  Although  the  pulse  amplitude  and 
intervals  of  transmission  are  fixed,  multipath  returns  bring  randomness  to  non-intended  receivers, 
including  GPS.  On  the  other  band,  the  frequency  range  of  motor  ignition  noise  is  approximately  100  MHz 
-  10  GHz  [19,20],  which  overlaps  with  the  GPS  band. 
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One  of  the  main  contributions  of  this  chapter  is  to  show  the  effect  of  non-Gaussian  noise  on  the 
GPS  receiver  performance,  specifically,  its  delay  lock  loop  (DLL).  We  consider  impulsive  noise, 
described  by  three  different  canonical  models,  namely,  Middleton  noise  [21],  generalized  Cauchy  noise 
[22]  and  symmetric  a-stable  (SaS)  noise  [23],  We  also  examine  the  automotive  ignition  noise  and  UWB 
noise  models  based  on  the  empirical  data  provided  in  [24,25]. 

It  is  well  known  that  DLL  is  the  main  component  of  the  GPS  receiver  which  performs  acquisition 
and  provides  the  pseudo  ranging.  The  DLL  correlates  the  incoming  digital  intermediate  frequency  (IF) 
signal  with  the  C/A  reference  code  from  three  different  shifting  positions,  referred  to  as  early,  late  and 
punctual.  The  results  are  then  sent  into  a  discriminator  for  arithmetic  processing.  Commonly,  the  system 
asserts  the  synchronization  when  the  discriminator  output  is  null.  However,  the  correlator  noise 
components,  which  result  from  the  correlation  between  the  incoming  noise  sources  and  the  early,  late  and 
punctual  C/A  reference  codes,  lead  to  a  tracking  error  and  loss  of  accurate  synchronization.  In  this 
chapter,  we  evaluate  the  impact  of  the  non-Gaussian  noise  on  the  DLL  discriminator  performance, 
incorporating  the  effect  of  the  front-end  bandpass  precorrelation  filter  and  sampling  rate.  The 
precorrelation  filter  reduces  the  out-of-band  noise  and  interference,  but  also  increases  the  correlation 
between  noise  samples,  which  in  turn,  increases  the  discriminator  tracking  error.  Both  the  tracking  error 
variance  and  the  discriminator  error  variance  are  derived  and  shown  to  be  highly  dependent  on  the  signal- 
to-noise  ratio  (SNR),  sampling  rate,  and  precorrelation  filter  bandwidth.  We  utilize  the  fact  that  due  to  the 
long  time-averaging  typically  performed  in  GPS  receivers  using  a  large  number  of  C/A  chips,  the  central 
limit  theorem  (CLT)  can  be  applied  to  characterize  the  correlator  output  noise  components,  which  are 
shown  to  be  Gaussian. 

Section  2  provides  an  overview  of  commonly  used  and  known  non-Gaussian  noise  models.  The 
precorrelation  filter  and  its  effect  on  the  DLL  discriminator  performance  are  discussed  in  Section  3. 
Requirements  on  the  filter  bandwidth  and  sampling  rate  are  also  established.  The  expressions  of 
discriminator  errors  and  their  statistical  characteristics  in  noise  are  derived  in  Section  4  in  terms  of 
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sampling  rate  and  precorrelation  filter  bandwidth.  Computer  simulations  are  provided  in  Section  5  and 
compared  with  the  theoretical  results. 


2.  Non-Gaussian  Noise  Models 


In  this  section,  we  briefly  discuss  several  commonly  used  non-Gaussian  noise  models.  We  present  the 
noise  probability  density  function  (PDF)  and  show  corresponding  sample  realization.  These  noise  sources 
are  employed  in  the  GPS  receiver  analysis  and  simulations  throughout  the  section. 

A.  SaS  Noise  Model 

The  a-stable  model  [23,26],  which  exhibits  sharp  spikes  and  occasional  bursts,  is  typically  used  for 
describing  the  impulsive  phenomena.  The  standard  SaS  PDF  is  given  by 

f(z)  =  ~Zj  .  ,  r(ak  +  l)|z|  sin[— — ]  ,  0<a<l  (1) 

7i  k=1  k!  2 


f  (z) = — jr  <JL  r(^±i)z2k 

na "  2k! 


1  <a<2 


The  parameter  a  assumes  values  in  the  range  (0,  2]  and  measures  the  “thickness”  of  the  tails  of  the 
distribution.  A  small  value  of  a  implies  considerable  probability  mass  in  the  tails  of  the  distribution, 
which  means  “severe  impulsive”.  There  are  two  important  special  cases;  the  PDF  is  Gaussian  for  a  =  2, 
and  Cauchy  for  a  =  1.  Figure  1(a)  shows  the  PDF  of  SaS  noise,  whereas  Figure  1(b)  presents  a  sample 
noise  sequence  for  a  =  0.5. 

B.  Middleton  Noise  Model 

Middleton  impulsive  noise  model  [21]  is  composed  of  both  Gaussian  and  impulsive  noise  components. 
The  probability  density  function  (PDF)  is  defined  by  two  parameters  A  and  r .  The  term  A  represents  the 
product  of  the  average  duration  of  the  pulse  and  the  average  number  of  pulses  occurring  in  unit  time. 
Small  values  of  A  increase  the  impulsiveness  of  the  noise,  whereas  large  values  of  A  move  the  model 

closer  to  a  Gaussian  distribution.  The  parameter,  r  =  aG2  /Oj2 ,  is  defined  as  the  Gaussian-to-impulsive 


44 


2  2  ‘ 

noise  power  ratio,  where  oG  represents  the  Gaussian  noise  power  and  a,  is  the  impulsive  noise  power. 

The  total  noise  power  is 


2  2,2 

o  =oG  +0t 


and  the  PDF  is  given  by 


p(z)=2 


e'AAm  1 


m=0 


m! 


:exp(--^— ) 


T,2  2°" 


(3) 


(4) 


with 


(m/A)  +  r 
i+r 


(5) 


If  A  is  sufficiently  small,  then  the  model  can  be  simplified  by  only  retaining  the  dominant  terms  in  (4) 
corresponding  to  m  =  0,  1,  2  [27].  In  this  case,  the  noise  PDF  function,  shown  in  Figure  2(a),  can  be 
approximated  by 


Pz(z)=  Z 


-  e‘AAm 


m=0 


m! 


^2na, 


r  exp(-  * 


2a. 


(6) 


The  above  PDF  is  an  even  function  with  heavy  tails,  as  illustrated  in  Figure  2(a).  A  sample  of  Middleton 
noise  sequence  is  depicted  in  Figure  2(b). 

C.  Generalized  Cauchy  Noise  Model 

A  generalized  Cauchy  probability  density  function  [22]  is  defined  in  terms  of  three  parameters,  a2,  k  >  0 
and  u>0. 


where 


p(z)  = 


B 

u  A 


(7) 


2  rq/k)  ,/2  B  ku-,/kr(n+i/k) 
r(3/k)J  ’  2AT(t))r(l/k) 


(8) 


In  the  above  equations,  T(  )  is  the  Gamma  function,  given  by 
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(9) 


T(a)  =  Jxa-Ie-Xdx 
0 

The  parameter  k  controls  the  impulsiveness  of  the  noise,  whereas  v  defines  the  noise  variance  and  o2  is  a 

scale  parameter.  In  the  example  shown  in  Figure  3(a),  we  set  k  =  0.2,  v=  40,  and  a  -  4.7.  Figure  3(b) 
shows  the  sample  sequence  of  the  noise. 

D.  Ignition  Noise  Model 

The  motor  ignition  noise  model  deduced  from  experiment  data  [24]  is  a  statistical  process  consisting  of 
two  distributions:  a  Weibull  distribution,  due  to  the  high  power  peaks,  and  a  Gaussian  distribution,  due  to 
low  power  values.  These  distributions  are  shown  in  Figure  4(a)  and  4(b),  respectively.  The  Weibull 
distribution  PDF  is  given  by 

p(z)  =  ab~aza~'e-(z/b)2  (10) 

We  set  a  =  1.14  and  b  =  4.00  in  Figure  4(a).  The  inter-arrival  times  between  successive  noise  peaks  are 
shown  to  be  exponentially  distributed,  as  illustrated  in  Figure  4(c).  A  sample  noise  sequence  is  shown 
in  Figure  4(d).  The  middle  strip  in  the  figure  corresponds  to  Gaussian  distribution  and  the  spikes 
are  weibull  distributed. 

E. UWB  Noise  Model 

The  UWB  signal  has  been  widely  surveyed.  The  signal  can  be  cast  as  impulsive  noise  since  it  consists  of 
intermittent  pulses  [28-30].  The  statistics  of  UWB  type  of  noise  generated  using  multiple  UWB  signal 
sources  are  reported  in  [25]  based  on  several  data  measurements.  With  the  assumption  that  the  phases  of 
the  band-limited  UWB  signals  are  uniformly  distributed,  it  was  concluded  that  the  amplitude  of  the 
aggregate  UWB  signals  approximates  a  Rayleigh  distribution,  which  is  shown  in  Figure  5(a).  Figure  5(b) 
gives  the  sample  sequence  of  the  noise  with  unit  variance. 

F.  Noise  Radar  Waveform  Model 

An  ultra-wideband  radar  is  a  radar  whose  fractional  bandwidth  relative  to  the  carrier  or  the  center 
frequency  is  greater  than  0.25.  It  transmits  ultra-wideband  noise  or  noise-like  waveform  towards  the 
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target  [12].  The  advantage  of  this  technique  lies  in  the  relative  immunity  of  the  random  noise-like 
unpolarized  transmit  signal  from  detection  and  jamming.  The  random  noise  is  modeled  as 

N(t)  =  A(t)cos[2jt(f  +  Af)tJ  (11) 

Typically,  A(t)  is  Rayleigh  distributed  and  Af  is  uniformly  distributed  between  [-0.5GHz,  +0.5GHz],  The 
center  frequency  of  the  noise  waveform  is  set  to  1.5  GHz  [12],  It  is  assumed  that  A(t)  and  Af  are 
uncorrelated  random  processes.  The  amplitude  PDF  of  the  noise  is  shown  in  Figure  6(a),  and  a  sample 
realization  is  shown  in  Figure  6(b). 

3.  Sampling  and  Filtering  Effects 

In  this  section,  we  discuss  the  effects  of  sampling  rate  and  filter  bandwidth  on  the  GPS  receiver. 
Conditions  on  both  variables  for  proper  tracking  and  synchronization  are  established. 

The  main  operation  for  code  synchronization  in  the  DLL  is  the  cross-correlations  performed 
between  the  received  data  and  the  receiver  reference  code.  Figure  7  shows  three  pairs  of  correlators 
required  to  produce  three  in-phase  components,  IE,  IP  ,  IL  and  three  quadrature  components  QE,  Qp, 
Ql,  respectively,  corresponding  to  the  early  (E),  punctual  (P)  and  late  (L)  reference  C/A  codes.  These 
codes  are  generated  from  three  shifting  code  positions,  but  with  the  same  intervals.  By  processing  the 
above  six  components,  the  receiver  may  construct  at  least  three  different  DLL  discriminators  [31], 
namely, 

Coherent 

D  =  (Ie  -IL)sign(IP)  (12) 

Early-minus-late  power  (noncoherent) 

D  =  (Ie2+Qe2)-(Il2+Ql2)  (13) 

Dot-product  ( noncoherent ) 

D  =  (Ie-Il)Ip+(Qe-Ql)Qp  (14) 
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where  sign(Ip)  is  the  sign  of  the  navigation  message  data  bit.  Synchronization  is  reached  by  locating  the 
correlation  peak.  Commonly,  the  discriminator  determines  the  peak  correlation  location  by  reaching  zero 
output  value.  The  output  of  correlator  j  (j  =  E,  P,  L),  corresponding  to  the  early,  punctual,  or  late  shifting 
positions,  can  be  written  as  [32]. 

Ij  =V2STR(Tj)cos<p  +  TiIj,  Qj  =V2STR(Tj)sin<p  +  T|Qj  (15) 

where  S  is  the  signal  power,  q>  is  the  residual  phase  error,  R(r.)  is  the  cross-correlation  function  between 
the  incoming  C/A  code  and  the  reference  code  corresponding  to  a  delay  x  -s  and  performed  over  T  (usually 
0.001)  seconds.  In  the  above  equation,  and  r)Qj  are  the  in-phase  and  quadrature  noise  components  of 
the  correlator  outputs.  Clearly,  if  g>  f  are  the  reference  code  sample  sequence,  and  n[jf  and  nQj, 
represent  the  input  noise  samples,  then  the  in-phase  and  quadrature  noise  components  can  be  expressed  as 

M  M 

Tlij=Zgj^nlj^  ’  1lQi=Zgj,^Qi^  (16) 

e=\  £=i 

where  M  =  K0T  is  the  number  of  samples  employed  in  the  computation  of  the  cross-correlation  (usually  it 
is  an  integer  multiple  of  1023),  and  K0  is  the  sampling  rate. 

A.  Sampling  Rate 

The  C/A  code  consists  of  1023  chips.  In  order  to  generate  the  proper  early-late  correlation  functions  for 
the  discriminator,  the  required  number  of  samples  per  chip,  K,  becomes  dependent  on  the  early-late 
correlator  spacing  d.  The  synchronization  position  typically  exists  at  the  midpoint  between  the  two  peaks 
of  the  early  and  late  correlation  functions.  Therefore,  in  order  to  accurately  find  the  synchronization  point, 
at  least  2  samples  should  be  generated  over  the  linear  range  of  the  early-late  correlator  spacing.  This 

requires,  at  minimum,  one  sample  for  each  half  correlator  spacing,  — .  If  we  express  the  spacing  as 

d  =  —  fraction  of  a  chip,  where  m  and  n  are  both  integers,  forming  an  irreducible  fraction,  then  half 
n 

spacing  is  — .  Therefore,  the  integer  number  of  samples  per  chip,  K,  should  be  no  less  than  2n,  if  m  is 
2n 
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not  a  multiple  of  2,  while  K  should  be  no  less  than  n,  if  m  is  a  multiple  of  2.  For  instance,  if  the  early-late 
correlator  spacing  is  3/5  chips,  then  the  number  of  samples  per  chip  should  be  at  least  10,  where  as  if  d  = 
2/5  chips,  then  m  =  2,  and  K  should  be  at  least  equal  to  5.  Details  of  this  example  are  shown  in  Figure 
8(a),  wherein  the  early-late  correlator  functions  with  spacing  of  2/5  chips  can  be  implemented  by  the 
sampling  rate  of  5  samples  per  chip.  The  discriminator  function  shown  in  Figure  8(b)  describes  the 
discriminator  output  due  to  the  tracking  error.  It  is  the  difference  between  the  early-late  correlation 
functions.  Without  multipath  and  noise,  the  discriminator  output,  assumes  zero  value  at  synchronization, 
which  can  be  estimated  by  two  reference  points  and  by  utilizing  the  linearity  over  the  delay  error  range  of 
the  discriminator  function.  In  Figure  9  (a)  and  (b),  we  reduce  the  sampling  rate  to  3  samples  per  chip.  It  is 
evident  that  in  this  case,  we  could  not  estimate  the  correct  synchronization  point,  since  the  sampling  rate 
doesn’t  provide  sufficient  resolution  for  holding  the  two  reference  points  in  the  linear  range. 

B.  Precorrelation  Filtering 

A  precorrelation  filter  is  applied  to  filter  out  noise  and  interference  outside  the  GPS  bandwidth.  Code 
tracking  accuracy  in  white  noise  has  been  examined  for  an  infinite  precorrelation  bandwidth  [31].  Filters 
with  finite  bandwidth  cause  correlation  changes  of  the  received  signals  and  remove  the  sidelobes  of  C/A 
code  spectrum.  Figure  10  shows  the  GPS  frequency  spectrum  due  to  the  application  of  a  Butterworth 
filter  of  2  MHz  bandwidth.  Typically,  the  filtering  adds  a  delay  to  the  incoming  C/A  code  relative  to  the 
local  C/A  code.  The  direct  result  is  the  offset  of  the  peak  location  of  the  correlation  function.  Figure  1 1 
shows  the  cross-correlation  function  between  the  filtered  C/A  code  and  the  local  original  C/A  code.  It  is 
observed  that  the  function’s  peak  location  does  not  correspond  to  zero  time  delay  error,  yielding  a 
tracking  error.  In  order  to  compensate  for  the  delay,  one  solution  is  to  subtract  the  known  induced  filter 
delay  from  the  overall  delay.  However,  even  with  delay  correction,  and  as  evident  from  Figure  11,  the 
correlation  function  due  to  filtering  is  not  absolutely  symmetric.  This  leads  to  a  slight  deflection  of  zero- 
delay  error  point  in  the  discriminator  function,  which  is  shown  in  Figure  12.  For  accuracy,  the  induced 
deflection  should  also  be  corrected. 
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One  way  to  mitigate  the  filter  effect  is  to  apply  the  same  precorrelation  filter  to  the  local  C/A 
code.  The  correlators  of  DLL  would  now  process  the  incoming  and  local  filtered  C/A  codes.  Figure  13 
shows  the  C/A  code  autocorrelation  function  after  applying  the  precorrelation  filtering  using  different 
values  of  the  filter  bandwidth.  The  autocorrelation  function  is  symmetric  with  no  delay  bias.  It  is  clear 
from  Figure  13  that  the  autocorrelation  function  changes,  depicting  a  more  smoothed  shape  for  narrower 
filter  bandwidth,  and  as  much,  compromising  the  linear  range  of  the  discriminator. 


4.  GPS  DLL  Performance  under  Different  Noise  Sources 


The  precorrelation  filter  also  modifies  the  correlation  properties  of  the  noise,  when  added  to  the  received 
signal.  The  ideal  band-limited  filtered  noise  has  a  sine  autocorrelation  function,  which  in  baseband,  is 
given  by 


where  Bf  is  the  precorrelation  filter  bandwidth.  In  Figure  14,  for  Bf  =lMHz,  the  noise  samples  spaced 


equal  to  or  larger  than  1  microsecond  (the  chip  width)  are  either  uncorrelated  or  weakly  correlated.  In 
essence,  oversampling  causes  the  sampling  interval  to  be  smaller  than  the  chip  width,  yielding  correlated 
noise  samples.  The  latter  adversely  affects  the  discriminator  performance  and  forces  the  discriminator 
output  to  assume  zero  values  at  the  wrong  signal  propagation  delay,  i.e.,  incorrect  satellite  range.  Figure 

15  shows  an  example  of  the  distorted  correlation  functions  and  corresponding  discriminator  function  of 
the  C/A  code  with  additive  filtered  generalized  Cauchy  noise  with  SNR  =  -10  dB. 

In  order  to  evaluate  the  discriminator  performance  in  noise,  the  variance  expressions  for  different 
noise  components  of  the  correlator  outputs  should  be  derived.  A  data  sample  in  the  noise  sequence  can  be 
considered  as  a  random  variable.  After  passing  through  the  precorrelation  filter,  the  noise  samples  change 
both  their  distributions  and  correlation  functions,  depending  on  the  filter  bandwidth.  For  instance,  Figure 

16  shows  the  PDF  and  the  sample  sequence  of  the  filtered  Middleton  noise.  Compared  with  the  original 
PDF  of  Figure  2,  it  is  clear  that  the  impulsiveness  and  density  have  been  reduced  due  to  filtering. 


50 


However,  the  uncorrelated  noise  samples,  whose  spacing  is  equal  to  the  half  of  the  main  lobe  of  the 
autocorrelation  function,  are  still  i.i.d.  (independent,  identical  distribution). 

To  derive  the  noise  error  variance  at  the  filter  output,  we  group  the  noise  samples  into  several 
sequences,  each  includes  uncorrelated  noise  samples  separated  by  an  inverse  of  the  filter  bandwidth.  The 
correlator  output,  therefore,  is  the  sum  of  the  correlations  between  the  receiver  reference  C/A  code  and 
the  different  grouped  noise  sequences.  It  is  noted  that  due  to  long  time  averaging  imposed  by  the 
correlation  loop,  the  Central  Limit  Theorem  (CLT)  can  be  applied,  yielding  a  Gaussian  random 
distribution  for  the  output  of  each  noise  sequence.  As  a  result,  for  the  sampling  rate  of  K  samples  per 
chip,  the  above  correlation  process  yields  K  Gaussian  random  variables,  the  sum  of  which  is  also 
Gaussian.  Accordingly,  the  in-phase  and  quadrature  noise  components,  r)(j  and  qQj,  in  (16)  are  both 

Gaussian.  The  histograms  of  those  two  components  for  the  case  of  filtered  Middleton  noise  using  2 
samples  per  chip  are  shown  in  Figure  17.  The  Gaussian  shape  is  evident. 

We  consider  the  early-minus-late  power  discriminator.  From  (13)  and  (15),  the  discriminator 

output  is 


D  =  (Ie2+Qe2)-(Il2+Ql2) 

=  (V2STR(TE)coscp  +  r|,E)2  +  (V2STR(te  )  sin  <p  +  t]QE  )2  (17) 

-  (V2ST R(tl ) cos cp  +  rijL )2  -  (V2ST R(tl )  sin  cp  +  qQL  )2 

where 


and 


te  =xp  +d/2  , 


tl  =  tp  -d/2 


(18) 


R(TE)  =  R(TP+d/2)  =  l-(TP+d/2),  R(tl)  =  R(tp  — d/2)  =  l  +  (xp  -d/2)  (19) 
assuming  the  prompt  time  delay  tp  is  bounded  into  the  discriminator  linear  working  range.  When  D  =  0, 
the  tracking  error  p  can  be  expressed  as 


P  —  Tp  ^D=0 


_ -tla.)cos<p  +  (iiQE  -qQL)sin<p] 


8ST(l-d/2) 


(20) 
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M  =  1023K ,  Bs=Kxl06  (27) 

The  derivation  of  the  above  equations  is  given  in  APPENDIX  A. 

All  terms  in  (21)  are  highly  dependent  on  SNR,  defined  as  (SNR  =  101og10(S/N0Bf )).  For 

relatively  high  SNR,  the  term  (Efri^2])2  ~(E[tiietiil])2  is  relatively  small,  compared  with 
2ST(l-d/2)2(E[qffi2]-E[T|IET|EL]).  Thus,  (21)  can  be  simplified  to 
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Varfo)  -  ggI£zl(j)!<lniE !lzgnB»iaJ)  =  lng!lzMMELj 

""7"7"  '"'7  4ST 


8S2T2(l-d/2)2 


(28) 


The  above  approximation  is  valid  for  SNR  above  -30  dB,  which  is  typically  the  case  for  the  GPS.  If  we 
incorporate  (22)-(2 6),  then  (28)  becomes 


Var(p)  =  i;[dM+  g  (M -i)(l  +d - i / K) -S‘n0,I:B r  /B,) 

+  2  g  (M-00-i/K)S^^_Y(M_ixl_|KdH/K)^,/B,)1 


biBf  /Bs 


K-Kd 


i7iBf  /Bs 


1=1 


i7tBf  /Bs 


(29) 


The  first  summation  in  (29) 


T(M-i)(l  +  d-i/K)^^> 

i=\  mBf  /Bs 

~TM(l  +  d-i/K)^/B>> 

1=i  ijrBf  /Bs 


(30) 


In  the  above  equation,  M-i  is  approximated  by  M,  since  the  upper  limit  of  the  sum,  K-Kd-1,  is  very  small 
compared  to  M.  We  rewrite  (30)  as 


K-Kd-1 


]T  M(l  +  d-i/K) 


sin(i7iBf  /Bs)_  1  .  ,„K  sin(i7cBf /Bs) 


i-1 


ia:Bf/Bs  Ts  i=1 


X  M(l  +  d  — i/K)- 


ijtBf  /Bs 


(31) 


where  Ts  is  the  sampling  interval,  and  can  be  expressed  in  terms  of  the  code  chip  width,  Tc ,  as 


T  =1^=1- 
s  K  M 


(32) 


K 

If  we  set  i  =  —  t ,  then  (31)  can  be  written  as 


i  K-Kd-1 

—  Y  M(l  +  d  — i/K) 

X  ti  iJtBf  /B 


sin(i7rBf  /Bs)t  _  M  sin[toBfK/(BsTc)] 


T  * 


(l  +  d-t/Tc)- 


t7rBfK/(BsTc) 


dt  (33) 


Using  the  relationships  in  (27), 
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M  (0-<*)Tc  . 

—  I  (1  +  d  - 1  /  Tc ) 

=“  r)T-(u.d-t/Tc)sinloiB'/<lfr‘>idt 

T  j)  IttR .  /no6  T  ^ 


sin[t7tBfK/(BsTc)] 

trtBfK/(BsTc) 


Tt 

KM  jO-d)Tc  ^ 
3-c 

=  A!MK 


+  d  -  t/T  ) 


t7iBf  /(10  Tc) 
sin[t7iBf  /(1Q6Tc)] 
tJtBf  /(106Tc) 


(34) 


dt 


where  Aj  is  a  constant,  independent  of  K.  Similarly,  the  other  two  terms  of  the  summation  in  (29)  can 
also  be  simplified  as  A2KM  and  A3KM ,  where  A2  and  A3  are  constants.  Thus,  (29)  can  be  written  as 


Var(p)  =  ^[d  +  (A,  +  A2  +  A3  )K] 


(35) 


We  substitute  (24)  into  (35),  with  Pn  =  N0BfT/M , 


Var(p)  = . [d  +  (A,  +  A2  +  A3  )K]  =  -^-[d  +  (A,  +  A2  +  A3  )K]  (36) 

Therefore,  the  tracking  error  variance  increases  for  increased  sampling  rate.  This  change  is  linear. 

Another  common  approach  to  evaluate  the  discriminator  performance  deals  with  the 
discriminator  error.  In  this  case,  the  error  value  is  the  discriminator  non-zero  output  at  synchronization. 
Compared  with  the  tracking  error,  which  is  expressed  in  time,  the  discriminator  error  is  analyzed  in 
amplitude.  It  is  straightforward  to  show  that  the  variance  of  discriminator  error  can  be  expressed  as 

Var(DlI[>=0)  =  16ST(l  — d/2)2(E[i|IE2]-E[r|IETiIL])  +  8(E2[TiIE2]-E2[riIEriIL])  (37) 

By  examining  both  error  variance  expressions,  we  find  that  the  reduction  of  correlator  spacing  decreases 
the  error  variance.  The  error  variance  can  also  be  lowered  by  decreasing  the  sampling  rate.  However,  as 
discussed  in  the  previous  section,  the  sampling  rate  itself  is  lower  bounded  with  a  value  decided  by  the 
correlator  spacing.  Figure  18  reveals  the  bound  on  the  tracking  error  due  to  minimum  sampling  rate  at 
different  correlator  spacing. 
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5.  Simulation  Results 


In  this  section,  we  present  the  results  obtained  from  20,000  Monte  Carlo  trials,  generated  based  on  the 
non-Gaussian  noise  models  described  in  section  2.  The  objective  is  to  estimate  the  effects  of  noise,  the 
precorrelation  filter  bandwidth,  and  sampling  rate  on  GPS  receiver  performance,  specifically  on  the  early- 
late  discriminator  tracking  error  variance.  We  compare  the  analytical  expressions  of  the  tracking  error 
variance  derived  in  Section  4  with  the  simulation  results.  The  first  set  of  simulation  demonstrates  the 
relationship  between  the  precorrelation  filter  bandwidth  and  both  the  tracking  error  variance  and  the 
discriminator  variance.  The  second  set  of  simulations  shows  the  discriminator  statistics  in  terms  of  the 
sampling  rate. 

First,  we  compute  the  values  of  discriminator  tracking  error  variance  over  a  range  of 
precorrelation  filter  bandwidths  between  2  and  10  MHz.  The  sampling  rate  is  fixed  at  30  samples  per 
chip,  and  the  early-late  correlator  spacing  is  set  to  0.4  chips.  All  non-Gaussian  noise  models  discussed  in 
Section  2  are  considered  and  have  given  almost  identical  results.  The  SNR  is  set  to  -30  dB.  Figure  19 
shows  that  the  tracking  error  variance  increases  as  the  precorrelation  filter  bandwidth  increases.  The 
simulations  follow  very  closely  the  analytical  results  of  equations  (28)  and  (37).  Increasing  the  filter 
bandwidth  beyond  that  of  the  GPS  reduces  the  noise  correlation,  but  at  the  same  time  allows  more  noise 
into  the  discriminator.  These  competing  forces  are  responsible  for  almost  the  constant  behavior  over  the 
range  4-6  MHz.  However,  the  figure  shows  that  a  large  precorrelation  filter  bandwidth  generally 
compromises  the  discriminator  performance,  irrespective  of  the  noise  statistics.  We  also  simulate  the 
discriminator  error  variance  versus  the  precorrelation  filter  bandwidth,  which  is  shown  in  Figure  20.  The 
result  is  consistent  with  that  of  the  tracking  error. 

The  second  set  of  simulations  considers  the  noise  models  of  Section  2.  Figure  21  shows  the 
different  values  of  discriminator  tracking  error  variance  with  a  range  of  sampling  rates  between  5  and  30 
samples  per  chip,  under  SNR  of  -30  dB.  The  precorrelation  filter  bandwidth  is  fixed  as  2  MHz  and  the 
early-late  correlator  spacing  is  set  to  0.4  chips.  We  find  that  the  tracking  error  variance  increases  as  the 
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sampling  rate  increases.  The  analytical  and  simulation  curves  in  Figure  21  are  very  close  and  both 
demonstrate  linear  behavior.  Figure  22  shows  the  relationship  between  the  discriminator  error  variance 
and  the  sampling  rate,  which  gives  the  same  result. 

Figure  23  provides  the  tracking  error  variance  with  2  MHz  precorrelation  filter  bandwidth  and  30 
samples  per  chip  over  a  range  of  SNR  between  -50  dB  and  10  dB.  It  is  clear  that  the  simulation  result 
tracking  error  variance  is  very  close  to  its  theoretical  counterpart  of  (28)  at  SNR  above  -30  dB. 

6.  Conclusions 

This  chapter  has  dealt  with  the  effect  of  non-Gaussian  noise  on  GPS  receiver  performance.  It  was  shown 
that  long  time-averaging  featured  in  the  GPS  DLL  is  responsible  for  producing  Gaussian  noise.  We 
considered  the  error  variance  of  the  early-late  discriminator  which  is  induced  by  the  correlator  noise 
components  of  impulsive  noise,  UWB  signals,  and  noise  radar  signals.  Both  expressions  of  the  tracking 
error  variance  (time)  and  the  discriminator  variance  (amplitude)  in  noise  are  derived.  The  results  show 
that  the  DLL  discriminator  statistics  are  highly  dependent  on  the  precorrelation  filter  bandwidth,  the 
sampling  rate,  the  early-late  correlator  spacing,  and  the  SNR.  The  precorrelation  filter  plays  a  key  role  of 
changing  the  correlation  between  noise  samples,  which  subsequently  affects  the  error  variance  of 
discriminator.  Higher  sampling  rates  increase  noise  sample  correlations,  which  leads  to  increased 
discriminator  tracking  errors.  Lowering  the  number  of  samples  per  chip  is,  however,  bounded  by  the 
early-late  correlator  spacing.  These  competing  requirements  as  well  as  the  tradeoff  between  noise 
correlation  and  noise  power,  as  a  result  of  increasing  filter  bandwidth,  are  analyzed  and  verified  by 
computer  simulations. 


APPENDIX  A 

First,  we  discuss  the  statistics  of  the  correlator  noise  components  r|IE,  T]^,  r|IP ,  rjQE ,  t]ql  ,  and  t]qp 
from  (16).  We  note  that, 
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M 

E[nIE2]  =  E[(£nWg(f))2] 


/=! 


M  MM 

=  E[£  (n(£)g(£))2  +  2n(l)g(l)£  n(£)g(£)  +  2n(2)g(2)£  n(Qg(0 

t= 1  1=2  1=3 

+  -  +  2n(M  -  l)g(M  -  l)n(M)g(M)] 


(Al) 


where  n(/)  represents  the  noise  sample,  g(£)  is  the  C/A  code  sample,  and  M  is  the  total  number  of 
samples.  Since  n(/;)  and  g(()  are  independent,  (Al)  becomes 


E[Tlffi2]  =  E[(£n(C)g(tO)2] 

i=l 

MM  M 

=  £  E[(n(/?)2  ]E[g(£) 2  ]  +  2£  E[n(l)n(£)]E[g(l)g(0)  +  2£  E[n(2)n(^)]E[g(2)g(£)]  (A2) 


f=l  1=2 

+  -  +  2E[n(M  -  l)n(M)]E[g(M  -  l)g(M)] 


t= 3 


For  simplification,  we  use 

Rf  (r)  =  E[g(a)g(a  +  r)]  =  l-r/K  for  r  =  0,  1,  2, K-  1;  Rf(r)  =  0  forr  =K,  K+  1, M  (A3) 

to  approximate  the  correlation  function  between  the  filtered  and  the  reference  C/A  code,  and 

sin(ntBf  /Bs) 


Rfn  (r)  =  E[n(a)n(a  +  r)]Pn 


nrBf  /Bs 


(A4) 


to  approximate  the  noise  correlation  function.  Thus,  (A2)  becomes 

E[tiffi2]  =  Pn [M  +  2£ (M  - r)(l - r / K) Sm(ntBf  ' Bs } 


] 


ncBf  /Bs 

where  the  total  number  of  samples  M  =  1023k.  Similarly, 

E[t,e2  ]  =  Ef^2  ]  =  E[nip2  ]  =  E[t|qe2  ]  =  E[tijL2  ]  =  E[T]jp2] 

Moreover, 

M  M 

EtWnJ  =  E[(£n(/)g(/))(£  n(*)g(*  +  Kd))] 


(A5) 


(A6) 
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M  M 

=  E[£  n(£)g(£)n(£)g(£  +  Kd)]  +  E[n(l)g(l)£  n(f )g(f  +  Kd) 

e=i  e=2 

M 

+  n(2)g(2)£  n(*)g(*  +  Kd)  +  •  •  •  +  n(M  -  l)g(M  -  l)n(M)g(M  +  Kd)] 

£~3 

2  M-l 

+  E[n(2)g(2)n(l)g(l  +  Kd)  +  n(3)g(3)£  n(/)g(/  +  Kd)  +  •••  +  n(M)g(M)£  n(f)g(f  +  Kd)] 

fc=i  1 


M  M 

=£  E[n(/)ii(/)]E[g(Og(*  +  Kd)]  +£  E[n(l)n(()]E[g(l)g(f  +  Kd)] 

£=2 

M 

+  £  E[n(2)n(/)]E[g(2)g(*  +  Kd)]  +  •  •  •  +  E[n(M  -  l)n(M)]g[(M  -  l)g(M  +  Kd)] 

2  M-l 

+  E[n(2)n(l)]E[g(2)g(l  +  Kd)]  +  £  E[n(3)n(f  )]E[g(3)g(f  +  Kd)]  +  •  -  +  £  E[n(M)n(f  )]E[g(M)g(f  +  Kd)] 

i=l  f=l 


K-Kd-l 


=  Pn{(l-d)M+  £  (M  —  r)(l - d - r / K) 


r=l 


sin(ntBf  /Bs) 
raBf  /Bs 


K1S_1  ,  .  sin(niBf  /B.) 

+  £  (M  -  r)(l  -  |Kd  -  r|  /  K) 


r=l 


} 


(A7) 


where  d  is  the  early-late  correlator  spacing  measured  in  chips. 

Next,  we  derive  the  expression  of  discriminator  tracking  error  variance,  assuming  the  tracking 
error  is  confined  into  the  linear  working  range  of  discriminator  function.  From  (15), 


D  =  (Ie2+Qe2)-(Il2+Ql2) 

=  (V  2STR(x  e  )  cos  (p  +  T|  [E ) 2  +(-v/2STR(xE)sin(p  +  T]QE)2 

- (V2STR(tl)cos<p  +  tJjl)2  -(V2STR(xL)sin<p  +  T|QL)2  (A8) 

=  {V2ST  cos((p)[l  -  (p  +  d  /  2)]  +  r\m  }2  +  {V2ST  sin(<p)[l  -  (p  +  d/  2)]  +  i)QE  }2 
- {V2ST cos((p)[l  +  (p-d/2)]  +  t]IL}2  - {V2Sfsin(9)[l  +(p-d/2)]-MlQL}2 

where  p  is  the  time  delay  of  C/A  code  autocorrelation  function  and  d  is  the  early-late  correlator  spacing. 

When  D  =  0,  the  tracking  error 

_  +t1qe  —t1  il  ’"‘Hql  +  2V2ST(l“d/2)[(r|J[E  —  t|il)cos(p  +  (t|qE  -“T|QL)sin(p] 
p_  _____  - 

Because  of  the  zero  mean  of  all  the  correlator  noise  components  and  (A6), 

E[p]  =  0 


(A9) 


(A10) 
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Thus,  the  tracking  error  variance 


Var(p)  =  E[p2]-(E[p])2  =E[p2]  (All) 

We  have  validated  that  all  correlator  noise  components  r|1E,  t|IL,  r)QE  and  r|QL  are  Gaussian.  Taking 
advantage  of  Gaussian  Joint-Variable  Theorem, 


E[Tllll21l31l4]  =  E[1llTl2]E[Tl3rl4]  +E[’ll^3]E[Tl21l4]  +  E[TllTl4]E[1l21l3] 
Using  the  above  expression,  it  is  straightforward  to  show  that 

Var(p)  =  _  ~ +2ST(l-d/2)2(E[r|IE  l-Eft^T]^]) 

8S2T2(l-d/2)2 


(A  12) 


(A  13) 


Finally,  we  present  the  derivation  of  the  discriminator  error  variance  expression.  At  the  point  of 
synchronization,  we  have 


R(xE)  =  R(xL)  =  l-d/2;  (A14) 

From  (15),  the  discriminator  output  is 

D1p=0=(Ie2+Qe2)-(Il2+Ql2) 

=  (V2STR(te  )  cos  9  +  t!1E  )2  +  (>/2STR(xe  )  sin  9  +  1]QE  )2 

-  (a/ 2STR(xl  )  cos  9  +  t|il ) 2  -  (V2STR(xL )  sin  9  +  r)OL  )2 

, -  ,  - _  (A  15) 

=  [V2ST  cos  9(1  -  d  /  2)  +  rllE  ]2  +  (V2ST  sin  9(1  -  d  /  2)  +  iiQE  ]2 

-  [V2ST  cos 9(1  -  d / 2)  + 1]^]2  - [V2ST  sin  9(1  -  d /  2)  +  riQL  ]2 

=  11ie2  +t1qe2  -t1il2  _t1ql2  +2-v/2ST(l-d/2)[cos9(TiIE  -%,)+sin<p(tiQE  —  tjql)] 

It’s  evident  that 


E[Dlp=0]  =  0  (A  16) 

The  discriminator  error  variance  is 

Var(p)  =  E[D  lp=0  ]  -  (E[D  lp=0])2  =  E[D  lp=02  ] 

=  16ST(l-d/2)2(E[t|IE2]-E[T]IET]IL])  +  8(E2[T]IE2]-E2[riIET]1L3) 
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Figure  3  (a)  PDF  of  generalized  Cauchy  noise  model,  (b)  Sample  sequence  of  Generalized  Cauchy  noise  with  unit  variance 
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Figure  4  (a)  PDF  of  the  power  peaks,  (b)  PDF  of  the  low  power  values, 

(c)  PDF  of  the  inter-arrival  times  between  successive  peaks,  (d)  Sample  sequence  of  ignition  noise  with  unit  variance 
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Figure  5  (a)  PDF  of  the  amplitude  of  aggregate  UWB  signals,  (b)  Sample  sequence  of  UWB  noise  with  unit  variance 
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Figure  6  (a)  PDF  of  the  amplitude  of  noise  radar  signals,  (b)  Sample  sequence  of  radar  noise  with  unit  variance 
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Figure  7  GPS  DLL  cross-correlation  process 
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Figure  10  The  frequency  spectrum  of  C/A  code  with  2  MHz  bandwidth  Butterworth  precorrelation  filtering 

(a)  before  filtered,  (b)  after  filtered 
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Time  delay  error  (chips) 

Figure  1 1  Cross-correlation  between  filtered  C/A  code  and  original  C/A  code 


Figure  12  (a)  Early-late  correlation  functions,  (b)  Discriminator  function 
with  1  chip  spacing  and  2  MHz  precorrelation  filtering 


Figure  13  C/A  code  autocorrelation  function  after  precorrelation  filtering 
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Figure  14  Correlation  between  noise  samples  with  1  MHz  precorrelation  filtering 


Figure  15  (a)  Early-late  correlation  functions,  (b)  Discriminator  function 
with  1  chip  spacing  and  2  samples  per  chip  under  -10  dB  impulsive  noise 


Figure  16  (a)  PDF  of  the  filtered  Middleton  noise,  (b)  Sample  sequence  the  filtered  Middleton  noise 
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Figure  17  (a)  PDF  of  correlator  in-phase  noise  samples,  (b)  PDF  of  correlator  quadrature  noise  samples 


Figure  18  Minimal  tracking  error  with  2  MHz  bandwidth  precorrelation  filter 


Figure  19  Discriminator  tracking  error  variance  through  different  bandwidth  precorrelation  filter 

under  non-Gaussian  noise 
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Figure  20  Discriminator  error  variance  through  different  bandwidth  precorrelation  filter 
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Figure  21  Discriminator  tracking  error  variance  with  different  sampling  rate 
under  non-Gaussian  noise 
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Figure  22  The  discriminator  error  variance  with  different  sampling  rate 
under  non-Gaussian  noise 


70 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188), 
1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any 
penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  0MB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

15/09/2005  Final 

3.  DATES  COVERED  (From  -  To) 

June  2004  -  September  2005 

4.  TITLE  AND  SUBTITLE 

Array  Processing  for  Interference  Suppression  in  GPS  Receivers 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

N00014-04-1-0607 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Amin,  Moeness  G.  (PI) 

i 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Villanova  University  -  Center  for  Advanced  Communications 

800  Lancaster  Avenue  -  Tolentine  1 19 

Villanova,  PA  19085 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

5-27725 

9,  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research 

875  N.  Randolph  Street 

1 2th  Floor,  One  Liberty  Center 

Code  02 

Arlington.  VA  222 1 7-5660 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

ONR 

11.  SPONSOR/MONITOR’S  REPORT 

NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  Public  Release;  Distribution  is  Unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

We  examine  the  performance  of  the  delay  lock  loops  (DLL)  in  GPS  receivers  in  the  presence  of  non-Gaussian  noise.  Interference 
and  noise  sources  for  GPS  receivers  may  assume  Gaussian  or  non-Gaussian  distributions,  specifically  when  operating  inside 
buildings.  Non-Gaussian  noise  may  contaminate  the  GPS  satellite  signals  and  compromise  the  receiver  delay  lock  loops  (DLL), 
producing  significant  tracking  errors.  These  sources  include  impulsive  noise,  ultra-wideband  (UWB)  signals,  and  noise  radar  signals 
for  target  tracking  and  through-the-wall  imaging  applications.  We  have  considered  non-Gaussian  noise  and  examined  its  effect  on 
GPS  correlation  and  discriminator  outputs  for  the  commercial  GPS  using  Coarse  Acquisition  (C/A)  code.  The  discriminator  hacking 
error  variance  is  derived,  incorporating  the  effect  of  noise,  the  front-end  precorrelation  filter,  and  the  sampling  rate.  Further, 
Interference  suppression  and  multipath  mitigation  in  Global  Navigation  Satellite  Systems  (GNSSs)  is  examined.  A  self-coherence 
anti-jamming  scheme  is  introduced  which  relies  on  the  unique  structure  of  the  coarse/acquisition  (C/A)  code  of  the  satellite  signals. 
15.  SUBJECT  TERMS 

Anti-jam  GPS;  Interference  Suppression;  Multipath  Mitigation;  Impulsive  Noise;  Antenna  Array  Processing;  Self-Coherence 


1  16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

OF 

PAGES 

Dr.  Moeness  G.  Amin 

u 

u 

u 

uu 

19b.  TELEPHONE  NUMBER  (Include  area  code ) 

610-519-4599 

Standard  Form  298  (Rev.  8/98) 


Prescribed  by  ANSI  Std.  Z39.18 


